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SONIC-BOX METHOD EMPLOYING LOCAL MACH NUMBER 
FOR OSCILLATING WINGS WITH THICKNESS 

By S. Y. Ruo 

Lockheed- Georg la Company 


SUMMARY 


A computer program has been developed to account approximately for the 
effects of finite wing thickness in the transonic potential flow over an 
oscillating Wing of finite span. The program is based on the original sonic- 
box program of Rodemich and Andrew and has been extended to include the 
effects of the swept trailing edge and the thickness of the wing. Account 
for the non-uniform flow caused by finite thickness is made by application 
of the local linearization concept. The thickness effect, expressed in terms 
of the local Mach number, is included in the basic solution to replace the 
coordinate transformation method used in the earlier work. Calculations were 
made for a delta wing and a rectangular wing performing plunge and pitch 
oscillations, and the results were compared with those obtained from other 
methods. An input guide and a complete listing of the computer code are 
presented. 


INTRODUCTION 


In reference 1, the sonic-box method computer program was developed for 
calculation of unsteady transonic flow aerodynamics for oscillating planar 
wings with unswept trailing edge by approximating the wing planform with a 
matrix of square boxes. Later, it was extended to include the swept trailing 
edge and control surfaces In reference 2. The sonic-box method uses a doublet 
velocity potential as the baste solution to satisfy the linearized transonic 
flow, unsteady small-perturbation velocity-potential equation with the asso- 
ciated boundary conditions. 

In references 3 and k, the wing thickness effect is partially recovered 
by the inclusion of local Mach number in the governing equation for the 



unsteady transonic flow. It uses the concept of local 1 inearlzation- to reduce 
the nonlinear smal 1 -perturbat ion equation to a linear one with non-constant 
coefficients. This Is further reduced to a linear equation with constant co- 
efficients by an appropriate coordinate transformation. This final equation 
and the associated boundary condition in the transformed space become identi- 
cal to those treated in the physical space by Rodemich and Andrew in reference 
1. The numerical results for the wing with thickness were obtained by adopt- 
ing the sonic-box method in the transformed space. Because of the assumptions 
made in deducing the governing equation to a manageable form, this technique 
is applicable only to relatively thin wings. That is, the local mean Mach 
number on the wing surface must not be very different from unity. Further, 
it is assumed that there is no flow separation and no strong shock waves on 
the wing surface. 

The computer programs developed In references 1, 2, and 3 use the least- 
square method to fit some of the input data, such as wing deflection or steady 
Mach number distribution on the wing, and to fit the computed velocity poten- 
tial with a form of pre-determined polynomial surface for the subsequent cal- 
culation of the unsteady pressure and the generalized aerodynamic force 
coefficients. The computer program described In reference 4 adopts the natural 
cubic spline for fitting calculated velocity potential and the spline-surface 
for fitting input modal deflections and Mach number distribution instead of 
the polynomial-surface fitting used in references 1, 2, and 3- The codes in 
references 3 and 4 allow the computation of generalized aerodynamic force 
coefficients for wings of zero and finite thickness; the swept trailing edges 
are allowed but not the control surfaces. 

The computer program presented in this report is developed according to 
the "alternate technique" described in reference 5. The coordinate transfor- 
mation technique as used in references 3- and, 4 fails when the, mean local, Mach 
number on the wing becomes very different from unity. In order to avoid this 
problem, an alternate technique was proposed in reference 5 to approximately 
account for the thickness effect, expressed in terms of mean local Mach number 
on the wing, by including It directly in the basic doublet solution to replace 
the coordinate transformation. The computer .program thus developed is smaller 
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than that of reference 4 and the amount of computation required has also been 
reduced. 

In general, the basic assumptions and limitations applied to the computer 
code in reference 4 also apply to the present code. However, the present 
formulation avoids the difficulties associated with the artificial wake and 
wing-surface fold-over due to multivalued transformation which limits, the 
usefulness of the coordinate transformation formulation of reference 4. The 
zero thickness wing portion of the computer code is unchanged from that of 
reference 4. Input is identical in both codes and the output differs very 
little between them. 


SYMBOLS 


b reference length (dimension = L) 

Cp pressure coefficient 

exp, e exponential function 
i ^ 

k reduced frequency, wb/Uoo 

L unit of length 

L|j generalized aerodynamic force coefficient 

M local Mach number 

T unit of time 

Ueo reference velocity { freest ream) , (dimension - L/T) 

x,y,z dimensionless Cartesian coordinates (reference length = b) 

6ij phase angle of Lij 

T maximum thickness to root chord ratio 

T*q magnitude of oscillatory dimensionless small perturbation velocity 

potential 

a angular velocity (dimension = radian/T) 

( subscripts denote quantity at leading edge 

( )|.g subscripts denote quantity at trailing edge 


3 



METHOD 


The computer program described in reference >4 is based on. the coordinate 
transformation technique to reduce the locally linearized equation with non- 
constant coefficients for nonzero thickness wing at sonic speed to a lin'ear 
one with constant coefficients. This linear equation with the associated 
boundary conditions can, then, be solved with sonic-box method. When the mean 
local Mach number on the wing becomes very different from that of the free- 
stream, the transformation may become multivalued and consequently an arti- 
ficial wake or wing-surface fold-over may be created in the transformed space. 
This technique fails once it happens. In order to avoid this problem, an 
alternate technique was proposed in reference 5 to approximately account for 
the thickness effect, expressed in terms of mean local Mach number on the wing, 
by including it in the basic doublet solution. 

The governing equation for unsteady transonic small perturbation velocity 
potential is 


'Poyy + - M^(2ikcPo^ - = 0 , 


( 1 ) 


where 


'Po(x,y,z) = T(x,y,z,t) *e , 
which is also equation (l) of reference 4. 

The basic solution for equation (l), representing a po.int doublet or iented 
parallel to the z-axis at the origin and satisfying the required condition at 
inf-lTi ity for a small finite region on the wing where the val.ue of M, the Mach 
number, is considered to be constant, may be written as 


0 , X <0 


fPo = 


ik zM^ 
2it 


exp {-iik[x + 


M2(y2+z2) 


]} , x>o, 


(2) 



in which M is regarded as a parameter. '.This solution satisfies equation (1) • 
only in .a small finite region of the wing; so. the solution may be considered 
to be of the locally linearized form. 

The only quantity in the program of reference A' requiring modification 
is the velocity influence coefficient for the wings with thickness. It is 
presently written as 

* - ^ 11 ^ - X [ ^ (3) 

E 

= ^ f exp [- -~{u + •^}] dudv {^) 

where H = length of the box side 

k = reduced frequency 
H = mean local Mach number - 

i = /:T • 

£ = kH 

E = box at (?,ri) 
u = (x-C)/H 
V = M(y~ri)/H 

« i 

The value of the velocity influence coefficient computed in the sonic-box 
computer program is with M = 1.0 in equations (3) and (4). Under this con- 
dition, the velocity influence coefficient is function of the wing geometry 
only. For M 1.0, it becomes function of the Mach number also. The value 
of the modified velocity influence coefficient required in this alternate 
technique to account for the wing thickness effect may be evaluated from the 
table computed for M = 1.0 condition for the same reduced frequency. 

To evaluate the modified velocity influence coefficient for this alter- 
nate technique, one may do the following: 

1. take the average value of the mean local Mach number at the center of the' 
receiving, (x,y) , and sending, (C,n), boxes, 
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2. multiply the spanwtse distance between these two box centers by the value 
of the average Mach number, 

3* interpolate the modified velocity influence coefficient from the original 
table for M = 1.0 with the value of the modified spanwise distance, v, 

4. multiply this value by the mean local Mach number at the center of the 
sending box. 

The rest of the computation remains practically unchanged except that the com- 
putation in the transformed space is totally eliminated. 

COMMENTS ON THE PROGRAM 

The velocity potential influence coefficients for a wing of zero thick- 
ness at a given frequency are only a function of the geometry. However, in 
addition to the geometry, they are also a function of the local Mach number 
distribution for the nonzero thickness wing under present formulation. It 
may be possible to perform the’ integration in equation (3) analytically with 
a new formula or with that already in the earlier program with some 
approximation. No attempt was made to derive- the totally new formulation. 

One of the approximate methods which was studied but not implemented in the 
present program is to ‘ subst i tute the local Mach number, M, in the integrand 
of equation (3) with (l-e), where e is a positive or a negative small number. 
After expanding the expotential function involving e term and neglecting 
all or higher terms, one obtains an approximate form of the integrand, 

for a doublet at the origin and z = 0, as follows: 

+ iek • exp (x + ^)] . (5) . 

I 

The exponential function in equation (5) is the same as that used In the case 
for M = 1.0 and the routines in the earlier sonic-box computer program may be 
utilized to perform the integration. Due to its complexity, and the addition- 
al ’computer storage and time required, this approximate method was not adopted 
to generate the new velocity influence coefficient matrix with the Mach ' • 
number effect. Instead, it is interpolated from the velocity influence 
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coefficient matrix for the zero thickness wing as described in the pre- 
ceding section. The Mach number appearing in the exponential function in 
equation (3) Is only associated with the distance between receiving and 
sending points and it is regarded as to modify the effective distance be- 
tween these -two points. Therefore, the average Mach number is used to 
maintain its interchangeability. Another Mach number in equation (3) is 
regarded as to modify the doublet strength. Since the integration is per- 
formed over the surface of the sending box, it is logical to use the Mach 
number at that point. This simplification in coupling the Mach number 
effect enables a reduction of the size of the computer code and the compu- 
tation time. The computed results appeared to be reasonable under the 
assumptions of small perturbation theory and local linearization concepts. 

In the present formulation, it implies as in the coordinate transforma- 
tion formulation (ref. A) that the Mach number variation in the spanwise 
direction is not large. The accuracy of these methods decreases when a large 
Mach number variation in spanwise direction exists. The present method, 
however, does not fail abruptly as does the coordinate transformation method 
when spanwise variation of Mach number becomes large enough to cause multi- 
valued .transformation and hence fold-over of wing-surface. 

Since no smoothing has been applied on either the input data or any com- 
puted values in the data fitting process during the computation, the calcu- 
lated unsteady pressure coefficient distribution may not be smooth and should 
be used with caution. In order to use It, the computed pressure coefficient 
should be put through a smoothing process such as the smoothing portion of 
the two-dimensional cubic-spline fitting routines in the present program. 

The pressure coefficient is obtained by differentiation of a set of numerical 
values whereas the generalized aerodynamic force coefficient is obtained by 
integration. Since integration itself is a smoothing process, the resulting 
generalized aerodynamic force coefficient is considered to be acceptable with- 
in the bounds of the accuracy of the numerical, techniques and the adequacy of 
the sonic-box method. The option of data smoothing is not provided in the 
three-dimensional spline-surface fitting process used in the present code for 
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input data such as wing defl ect Ions- and mean local Mach number. The spline 
surface Is required to go through all input points. 


RESULTS 


Sample calculations are made for a delta wing and a rectangular wing 
oscillating in plunge (Mode 1) and in pitch about the apex (Mode 2). The 
mean angle of attack is zero and the freestream is at sonic speed. 

Delta Wing 

The delta wing considered here is a flattened elliptic cone of aspect 
ratio 1.5 and thickness-to-root-chord ratio x = 0.1. Convergence with re- 
spect to the number of boxes along the root chord for the generalized force 
coefficients (L}j) due to plunge and pitch about the apex, at a reduced 
frequency of k = 0.2, is shown in figure 2. The maximum numerical difference 
within the applied range of 15 and 30 boxes along the root chord is about ^ 
percent, and the trend of convergence with and without thickness is essential- 
ly the same. Based on the results shown In figure 2, it appears that the gain 
in convergence by using a large number of boxes to represent the wing is not. 
obvious as compared with a fortuitous selection of the number of boxes to 
use. The numerical fluctuation in the convergence plot is largely caused by 
the box arrangement along the wing leading edge which, in turn, is dependent 
on the number of boxes selected for use along the wing root chord. Contri- 
bution from the partial boxes along the leading edge has been taken into 
account, but the fluctuation still exists. 

The variation of each force coefficient, using 30 boxes along the root 
chord, versus the reduced frequency is plotted in figure 3. The results from 
figure 7 of reference h are also shown. The numerical difference between the 
results for wings with and without thickness is very small, generally less 
than one percent. This is a result of .the Mach number, at each box-center 
used in the computation, lying within the narrow range of 0.92 and 0.98 in 
chordwise direction and remaining constant in spanwise direction (see fig. 8 
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of ref. A)., However, the thickness effect on flutter speed can be signifi-' 
cant (ref. 6). The results for the case with thickness obtained, from the - 
present method and that of reference, 4 are, not very different. 


Rectangular Wing 

The rectangular wing considered here has aspect ratio 2.0 and a biconvex 
(circular arc) airfoil with th ickness-to-chord ratio t = 0.0521. The vari- 
ation of each force coefficient, using 20 boxes along the root chord, versus 
the reduced frequency is plotted in figure h. The results obtained from the 
present method, and from references 4 and 7, are included in the figure. The- 
thickness effect on the rectangular wing is seen to be slightly larger than . 
that on the delta wing. This is probably caused by the wide range of Mach 
number variation (fig. 6) on the rectangular wing, even though the thickness 
ratio is only 0.0521 for the rectangular wing against 0.1 for the delta wing. 

The present method predicts values higher than either Landahl's results 
(ref. 7) for the zero thickness case or the results of reference 4 for the non- 
zero thickness case. The difference of the generalized aerodynamic force co- 
efficients for the nonzero- thickness case between the results obtained from 
the present method and that of reference 4 is quite large. This might be 
caused by the difference in interpretation of the effective distance between ■ 
the sending and the receiving points in the present method and the coordinate 
transformation method used in reference 4. It is felt, however, that the 
interpretation used in the present program is more physically sound than that 
used in reference 4. The phase angle predicted by the sonic-box method at 
very low reduced frequency becomes meaningless when the magnitude of any force 
coefficient approaches to zero with decreasing reduced frequency (for example, 
see figs. 4(b) arid 4(d)). This is due to numerical inaccuracy and not to any • 
inadequacy of the method. 

The steady-state pressure coefficient ^obtained from reference 8 for the 
rectangular wing considered here is shown In f-igure -5(a) for the chordwise 
(x-di recti on) distribution and in fig,ure 5(b) for the spanwise (y-d i rection) 
distribution. The corresponding Mach, numbers at the box-centers, interpolated,- 
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from the fitted spline-surface, are plotted in figures, 6(a) and 6(b)>. The.'? 
interpolated values deviate from the Input data .more near -the leading edge 
than near the trailing edge, this probabl.y was caused: by.-the use of more- 
dense spacing of input points near the leading edge as. .compared with those 
near the trailing edge in chordwise direction and by the lack of input 
points near the leading edge in spanwise direction, especially in the in- 
board portion of the wing. A better fit than that shown in figures 6(a) and 
6(b) may be obtained by using more evenly spaced, .input points, than those 
shown in figures 5(a) and 5(b). 


CONCLUDING REMARKS 


A sonic-box method computer program is presented for the application of 
a local 1 inearizat ion .concept capable of accounting approximately for wing 
thickness effects in unsteady sonic flow. The thickness effect, expressed in^ 
terms of the local steady Mach number, is directly included in the basic 
solution. The local doublet strength is adjusted from the sonic flow con- 
dition to that for the local flow, and the governing equation is reduced to 
the one used in the original sonic-box method for zero thickness wings. Thus, 
the original sonic-box method concept can be used directly to treat nonzero 
thickness wings. 

Convergence of the numerical results with respect to the number of boxes 
used in representing the wing planform seemed to depend more on the arrange- 
ment of the boxes along a swept leading edge than on the total number of 
boxes used, even though the partial .boxes along the leading edge were in- 
cluded in the computation. For a wing with unswept leading edge, the use of 
a small number of boxes (say, 15 to 20 along the root chord) appeared to be . 
sufficient to obtain results that were essentially converged. 

When the input data requi re- spl ine-surface -f I tt ing, the input points must 
be selected in such a way that they are as. uniformly spaced as possible to. 
avoid locally-concentrated' large errors. A smoothing option for the two- 
dimensional cubic-spline has been included In the present program, but it was 
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not utilized in the-sample runs shown in this report. Since the box method 
itself is numerical .in nature, the distribution of calculated values may not 
always be smooth; and it may become necessary to perform the smoothing before 
any gradients are evaluated. 


Based on the sample runs made, the contribution due to thickness was not 
found to be very large in comparison with the results calculated from the 
coordinate transformation method. Due to the lack of reliable experimental 
data, it is rather difficult to assess the validity of the present approach 
in accounting for thickness effects. 
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(c) Li ft due to pi tch. 


(d) Moment due to pitch. 
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APPENDIX 
COMPUTER PROGRAM 


The computer program listed In this appendix is dimensioned to handle 
a maximum of 30 boxes, either in chordwise or spanwise directions, in approxi- 
mating the wing planform. The maximum numbers of leading and trailing edge 
segments are, respectively, 7 and 2 per semispan. The program can handle up 
to 3 wing deflection mode shapes. The maximum number of points used in 
spline-surface fitting is 100, so .the maximum number of input points to 
describe the wing deflections and the Mach number or pressure coefficient 
distribution is also 100. These limitations can easily be increased by 
changing the dimensions of the corresponding variables in the computer program. 

in order to activate the smoothing option in two-dimensional cubic-spline 
data fitting, it is required that a two-digit fixed point number be assigned 
■to the last argument, NSMOS, of subroutine SPLNl . The right digit is for 
the control of the number of smoothings desired; and the left-digit is for 
pre-interpolation, zero for omitting pre-interpolation and non-zero for in- 
cluding pre-interpolation. The pre-interpolation is a process to increase the 
number' of known points to be used in the interpolation by inserting an addi- 
tional point betwen every successive pair of input points in the original set. 
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Input Guide 



Data are input through the'.subroutine DATRD using the. one d.imens lonai 
array DA with a size of 1005. The allowable maximum number for some of the 
input data as Indicated below may be changed if the dimensions of the corre- 
sponding storage array and computational operations are also changed accord- 
ingly. Subroutine DATRD initializes DA(l) through DA(22) to blank, the 

weighting factors in DA{104), DA (108), , DA(500) to 1.0, and the remaining 

portion of the DA array to 0.0.. Consequently, these are the default values. 
The layout of the array DA(k) as it is presently .used is similar to that in 
reference 4 and is as follows: 

1-7: Title 

8-12: Not used 

I 3 -I 9 : Mode title 

20-22; Not used 

23 : Frequency, (cycle/sec) 

2h: Overall length of wi.ng in streamwise direction, (ft or meter) 

25: Speed of sound of the freestream, (ft-/sec or meter/sec) 

26 : ( 0 ) - indicates the frequency is the first one for a new wing 
(l) - indicates the frequency is the additional one for the 
same wing 

27 : Number of boxes in streamwise direction (maximum 30) 

28 : Number of deflection modes (maximum 3) 

29: Number (m) of segments of leading edge per semispan to be 

given, excluding segment from origin to yo (m[^ 3 j 5 = 7) 

30-44: Coordinates of points on the’ leading edge, (ft or meter) 

(in sequence of yp .Xi ,Yi ,X 2 ,y 2 ,-:-,Xm,ym) . mmax = 7 
45 : Number (n) of segments of trailing edge per semispan to be 

given, (default: unswept trai 1 i ng, edge) , nj,^gj^ = 2 
46-48: Coordinates of points on the trailing edge, (ft or meter) 

(in sequence of xo,yi,xi for n = 2, 
or xq (only) for n = 1, 
no input for n = 0; 

last trailing edge point coincides with the last leading edge 
point and is set internally) 


25 





50 : 


51: 


52-53: 

54-70: 

71-72: 


73-95: 
96 : 


97: 


98 : 

99-100: 

101 - 500 : 

5 OI- 70 O: 

701 - 1000 : 


Number of boxes allowed for upstream Influence (if- this 
location is left blank or assigned a zero, it will assume 
DA(49)=DA(27) and in no case DA(49) >DA(27) is allowed). 

(0) - indicates to calculate cases with- and without thickness 

effect 

( 1 ) - indicates to calculate case without thickness effect only 

( 2 ) - indicates to calculate case with thickness effect only 
Indicator to suppress calculation of potential for a mode 

( 0 ) - no suppression 

(1) - suppression 

Coefficients of the deflection polynomial (in the sequence of 
ag and a^) 

Not used* 

Coefficients of the Mach number distribution polynomial (in 
the sequence of ag and ai) 

Not used* 

Indicator of the type of wing thickness effect input 

( 1 ) - pressure coefficient 

( 2 ) - Mach number 

Number of points at which pressure coefficient or Mach number 
to be given 

Number of points on which deflections to be given 
Not used* 

Deflection data for a maximum of 100 points (in the sequence 
of X, y, deflection and weighting factor) 

Not used** 

Pressure coefficient or Mach number data for a maximum of 100 
points (in the sequence of x, y and pressure coefficient or 
Mach number) 


The remaining part of DA array is used for the control of intermediate 
results print out. When the latter is desired, a non-zero positive integer 
number should be entered at locations in the DA array corresponding to the 
information from the one particular subroutine that is needed. 
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1003 : 


1002 : 


1003: 


1004: 

1005: 


original 

OF POOR quality 


CBA: for wing deflection (DRED) 

A=1, for NEW=1 


B=1 

C 

BA: 

A=1 

B 


for NEW=2 


spline-surface fitted results 


- not appi 1 cable 
for wing upwash (WVAL) 

- upwash 

- not appl icable 

FEDCBA for velocity potential (BOXP and BOXPO) 

A=1, for NEW=1 
B=1, for NEW=2 
C=l, for NEW=1 
0=1, for NEW=2 
E=l, for NEW=1 
F=1, for NEW=2 
A; for Mach number (MRED) 

A=1 - spline-surface fitted results 

DCBA: for wing shape (SHAPE and PLNFM) 

A=1 , for NEW=l) distributions of box, box area, leading and 
B=1 , for NEW=2j trailing edges 
C=1 - Mach number at box centers 

D - not appl icable 


velocity potential 

Influence coefficient and solution matrices 
pressure coefficient 


The format of the input data card is (A1 , A5, l6, 6A10, A8) . The first 
field is for the control of clearing the data- array, DA, for a new wing (+) 
and the control to indicate the end of the set of data (-) . The second field 
is the indicator for the type of data, either numeric (blank) or alphameric 
(ALPHA). The third field is the designator for the relative location in the 
data array of the first number to follow in the fourth field. If this field 
Is left blank, or a zero is entered, the execution will be terminated. The 
fourth and fifth fields are for five consecutive Input data each occupying 12 
columns plus 8 blank columns at the end. All the fixed point numbers are right- 
adjusted and the decimal point for the floating point number must be included, 
if an input datum is left blank, no change at the storage location for that 
particular datum in the data array will occur unless the set of the input 
data is for a new wing. 
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Those storages currently not used in array DA marked with * are reserved 
for future improvements in the method used for the functional form of data 
input,- Those marked wi.th ** are reserved, for the case, where large' numbers 
of data input points for deflections- or Mach numbers are required. 
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Sample input 


A typical input data deck set-up for an aspect ratio 1.5 delta wing 
having an elliptic lateral cross-section with 10^ thickness ratio performing 
plunge and pitch about its apex are given below. 

The input format is (A1 , A5, 16, 6A10, A8) . 

Card 1: title of the case under consideration. 

Card 2: title of the first mode of deflection. 

Card 3: first frequency (cycle/sec), center.l ine 'chord length (ft), 

reference velocity (ft/sec). 

Card k: number of boxes along the centerline chord, number of deflection 

modes, number of total leading edge segments of the wing. 

Card 5: spanwise coordinate (ft) of the first section of the leading edge, 

chordwise and spanwise coordinates (ft) of the next section (the 
sequence is yg, x^, y^ — e.g., see figure 1). 

Card 6; first trode of deflection f ~ 1.0 

the “ - " sign indicates the end of the group of data cards to 

be read at this stage. 

Card 7: title of the second mode. 

Card 8: second mode of deflection f = O.lx. 

Card 9: identification of the type of Input regarding the wing thickness 

effect (Mach number for this case), number of points on the wing 
this information to be given. 

Cards 10 to 69: 

chordwise and spanwise coordinates (ft) of a point on the wing, and 
the Mach number at this point. 

the " - " sign on the last card indicates the end of the group of 
data cards to be read at this stage. 

Cards 70 and 71 : 

additional frequencies for the same wing, one card is read in at 
one time. 

Card 72: blank card to make an exit from the computer. 
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The card Iroages are as follows: 


1 2-3 5 6 7 9 

1234S67890l23«567890la3a5579qol?3fc5fi79901Z345678901234S67890l23«567890l234S67890 


♦«LPHA 

laSPECT PATIO 1,5 

DELTA WING (TAU*0,10> 

1 

ALPHA 

13PLUNGE 



2 


23 0.159134941 

lo.fl 

1000.0 

3 


27 30 

2 

•1 

4 


30 0.0 

10. 0 

3.75 

5 

- 

52 1.0 



6 

ALPHA 

13P1TCH ABOUT root 

leading 

EDGE X»0,o 

7 


52 

0,0 

0.1 


8 

96 

1 

60 


9 

701 

0,1 


0,140232 

10 

704 

0.3 


0.13S125 

11 

.707 

0.5 


0.132175 

12 

710 

0.7 


0,129943 

13 

713 

0.9 


0.128084 

14 

716 

1.3 


0.124980 

15 

719 

1.7 


0,122340 

16 

722 

2.1 


0.119963 

IT 

725 

2.5 


0,117746 

18 

728 

2.9 


0.nS<>27 

19 

731 

3.3 


0.113564 

20 

734 

3.9 


O.llOSlO 

21 

737 

A.S 


0.107426 

22 

740 

5,1 


0.104243 

23 

743 

5.7 


0.109872 

24 

746 

6.3 


0.097215 

2S 

749 

6.9 


0,093129 

26 

752 

7.5 


0.088395' 

27 

755 

7.9 


0.084697 

28 

758 

8.3 


0,080342 

29 

761 

8.7 


0.074970 

30 

764 

9.1 


0,067825' 

31 

767 

9.3- 


0,063062 

32 

770 

9.5 


0.056823 

33 

773 

9.7 


0,047661 

34 

776 

9.9 


0.029651 

35 

779 

1.7 

0.5 

0,128340 

36 

782 

3.3 

1.0 

0,113564 

37 

785 

4.5 

1.0 

0.107428 

38 

788 

5,7 

1.0 

0.100872 

39 

7.91 

6.9 

1.0 

0,093129 

40 

7.94 

6.3 

r.'o 

O.OS0342 

41 

797 

9.5 

1.0 

0.056823 

42 

800 

4.5 

1.5 

0.107428 

43 

603 

6.3 

2.0 

0,097215. 

44 

806 • 

7,5, 

2.0 

0.088395 

45 

809 

8,7 

2.0 

0.074970 

46 

812 

9.7 

2.0 

0.047661 

47 

815 

8,3 

2.5 

0.080342 

48 

818 

9.7 

3.0 

0.047661 

49 

821 

2.5 

O.S 

0.117746 

50 

824 

3.9 

'0,5 

0.110510 

5) 

827 

5,1 

O.S 

0,104243 

52 

830 

6.3 

0.5 

0.097215 

53 

833 

7.5 

0.5 

0.088395 

54 

836 

9.1 

0.5 

0.067825 ' 

55 

839 

9.9 

O.S 

0.029651 

56 

842 

5.1 

1.5 

0.104243 

57 

845 

6.3' 

• 1.5 

0.097215 

58 

848 

7.9 

1.5 

0.004697 

59 

851 

9.1 

,1.5 

0.067825 . 

60 

854 

9.9 

1.5 

0.029651 

6] 

857 

5.7 

2.0 

0. 100872 

62 

860 

6.9 

2.5 

0.093129 

63 

863 

9.5 

2.5 

0.056823 

64 

866 

9.9 

2.5 

0.029651 

65 

869 

8.3 

3.0 

0.080342 ' 

66 

872 

9,1 

3.0 

0,067825 

67 

875 

9.5 

3,5 

0.056823 - 

68 

878 

9.9 

- 3;s 

0,029651 

69 

23 r. 

591549407 


1 

70 

23 6, 

366197628 


1 

71 

72 
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Sample Output 


OEIGINAL BAGJ) 1ft 
JWR QU^ 


ASPf«;T ■?«TIO 1.5 D:LTA WING ITAU=0.10I 

10 BOXES AEONS ROOT C10RQ ROOT CH1RD LENGTH 

REOOCEO FREQUENCY = .010 FREE' STREAN VELOCITY 

FREQUENCY = 1.592E-01 CYCLE/SEC 

HOOF NO. 1 PLUNGE 

HOOF NO. 2 PITCH ABOUT ROOT LEADING EDGE X=0.0 


GENF»aLIZFO forces (NC thickness EFFECT) 


POPFS 

ORES. OEFL. 
1 1 

1 2 


REAL PART 
2.153V5E-06 
-1.69639E-07 


IPA6 PART 
-2.92783E-02 
-1.61159E-02 


GENERALIZED FORCES <N1 THICKNESS EFFECT) 


"TOES 

ROES. CEFL. 
2 1 

2 2 


REAL PART 
-2.A2738F*00 
-1.61 162F»00 


IRAQ PART 
-2.AA689E-02 
-1.813n<>E-02 


»00P NO. 1 
"ODF NO. 2 


GENfOALIZEO FORCES IwITH THICKNESS EFFECT) 


•'I'rf s 

o=ES- DEFL. 
1 1 

1 2 


R'^AL PART 
-1 .00531E-05 
-1.29333E-05 


I«AG part 
-2.93535E-02 
-1 .61 53AE-02 


GENERALIZED FiRCES lUlTH THICKNESS EFFECT) 


** T-E S 

ORES. CEFL. 
2 1 

2 2 


REAL “ART 
-2.935«1E»00 
-1.61 539f*00 


IHAG PART 
-2.33020E-02 
-1.6889RE-02 


ABS. VALUE 
2.42783E-02 
1 .61 15OE-02 


A95. VALUE 
2.42800E+00 
1 .61173E+0') 


ABS. VALUE 
2.43535E-02 
1 .61534E-02 


ABS. VALUE 
2.43552E*00 
1 .61548E+00 


10-00 FT 
1000.00 FT/SEC 


PHASE ANGLE 
-69.9949 
-90.0006 


PHASE ANGLE 
-179.4226 
-179. 3553 


phase angle 

-90.0237 
. -90,0459 


PHASE ANGLE 

-179.4518 

- 179.4011 
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AS'’FCT SATIO 1.5 D5LTA UlMG (TAOrO.lOl 


JO ="XFS ALON5 ->T1T CHORD ROOT CHORD LFHGTH 

REDUCED FRFOUi^^Cy = .100 FREE STREAM VELOCITY 


FO fOUENCY 

= 1 

•59?f+00 CYCLE/SEC 

POOF 

NO. 

1 

“OOf 

NO. 

2 


oE"JERALIZED F'^RCES IN'- THICKNESS EFFECT) 


“OPES 

RRES. CEFL. 
1 1 

1 2 


REAL PART 
5.6A26PF-0A 
1.2C006E-0A 


IHAG PART 
-2 .A21A1E-01 
-1 .6060AE-01 


generalized FORCE.S (NC thickness EFFECT) 


“ODES 

P°ES. OEFL. 
2 1 

2 2 


REAL PART 
-2.«?556E»00 
-1.6094RF+00 


IMAG PART 
-2.A51 5JE-01 
-1.81827E-01 


ABS. VALUE 
2.A2I AlE-01 
1 .6060AE-01 


ABS. VALUE 
2- A3792E*00 
1 .61972E*00 


“OOE A'O. 1 
“ODE NO. 2 

generalized forces CUITH THICKNESS EFFECT) 


MODES 


PRES. 

OEFL. 

REAL PART 

I HAG PART 

ABS. VALUE 

1 - 

1 

-7.33272E-00 

-2 .02506E-01 

2.02507E-01 

1 

2 

-1 .03563E-03 

-1.6059PE-01 

1-.6 06 01E-01 


GENERALIZED FTRcES tUlTH 

THICKNESS EFFECT) 


MODES 

PRES. DEFL. 

REAL PART 

INAG PART 

ABS. VALUE 

2 

1 

-2.'n011F»00 

-2.30356E-01 

2.00138E+00 . 

2 

2 

-1.61039E»00 

-1.70175E-01 . 

1.61936E+00 


10.00 FT 
1000.00 FT/SEC 


PHASE ANSLF 
-89.9138 
-80.9572 


PHASE ANGLE 

-170.2287 

-173.5505 


PHASE ANGLE 
-90.1732 
-90,3695 


PHASE ANGLE 

-170.0915 

-173.9678 
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ORIGINAL PAG® i?. 
OF POOR QUALm 


AS'^fCT RATIO 1.5 OFLTA UIN6 <TAU=0.10) 



30 BOXES 

ALONG ROOT CHORD ROOT 

CHORD LENGTH = 

10.00 FT 


REDUCED FREQUENCY = , 

AOO FREE 

STREAH VELOCITY = 

1000.00 FT/SEC 


FREOUENCY 

= 6.366E*00 

CYCLE/SEC 




NODE 

NO. 1 





NODE 

NO. 2 





generalized FORCES CNC 

THICKNESS EFFECT) 



finES 

PRES. DEFU. 

REAL PART 

IHAG PART 

ABS. VALUE 

PHASE ANGLE 

I 

1 

2.60265E-02 

-9.51379E-01 

9.51735E-01 

-88.4330 

1 

2 

1.9370AE-02 

-6.28AA1E-01 

6.287A0E-01 

-88.2345 


GENERALIZED FORCES CND 

THICKNESS EFFECT) 



NODES 

PPES. DEFL. 

REAL PART 

IHAG PART 

ABS. VALUE 

PHASE ANGLE 

2 

1 

-2.A1 529E*00 

-9.RR983E-01 

2.61373E+00 

-157.5296 

2 

2 

-1.60030E»00 

-7.9A07AE-01 

1.76A83E*00 

-155.0635 


fOOF NO. 1 
«ODF NO. 2 

GENf'»ALIZFO FOfiCFS IWITH THICKNESS EFFECT) 

P"»rES 


PRES. 

DEFL. 

REAL PART 

INA6 PART 

ABS. VALUE 

PHASE ANGLE 

1 

1 

1.95356E-02 

-9.43362E-01 

9.43564E-01 

-88.8137 

1 

2 

I.17730E-02 

-6.19079E-01 

6.19191E-01 

-88.9105 ■ 


GENERALIZED FORCES (WITH 

THICKNESS EFFECT) 



NODES 

PRES. TEEL. 

REAL PART 

INAG PART 

ASS. VALUE 

PHASE ANGLE 

2 • 

1 

-2.3R926F*00 

-9.75339E-01 

2.58993E+00 

-157.3775 

2 

2 

-1.5?205E+00 

-7.16462E-01 

1.73672E*00 

-155.6357 
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Program Listing 




} . (l-HOt,jjTi'Ul,Tiv-t3 = lHPUT,TAPE6*=0UTPUTJ 

SB0XR2 

z 



nil - r *■ > 

SB0XR2 

3 

1 


JP iJxlS along THt CENTERLINE CHORD ALLOWED 

SB0XR2 

4 



STATtHtNT 

SB0XR2 

S 



-. --A>lr'..r sv.fiLi' jf nOLtS ALLOXEC IN DIMENSION STATEMENT 

$B0XR2 

6 

c 


.I’ .r .ijXtS allowed for SPLINE FITTING 

SB0XR2 

7 

L 


}lnc.r*jlOK STATEHENT 

S80XR2 

8 



.'/./ * / ( 3C 3 » t If { ju ) * xlE 1 JU 1 » xi E C 30 } 1 EOG ( 30) « AR i 30* 30) 

A 

2 



L ^r.Ll jU f 3-J ) ,"<lT(2»30 ) »MLC( 2» 30) .MB »AMA( 30*301 

A 

3 



.Xi.3MbJ,Yi:JG( -i) .NtSCI ,nS.XT0G(3) .YT0GI3) .NTDCI.NST 

A 

4 



C ,bFw>JU)0).srdy(lGG»>$F«H(103).KSFM.T(103.10AI.HB,NN 

A . 

5 



L .jAIICjS) ,A c.^A,CK,D,OH,oI,YhAX,JHAX,JMAK,IEOG.IH,L.HEH 

A 

6- 



L'lNiMSioN A(2, 60 1 jO ) . :> ( 2 1 30. 30 ) 

SBQXR2 

10 



. ,SFDi (IOC ,3) ,5 FDy“« 100, J).SFDH(103,3) .KSF0I3) ,G(3) 

SB0XR2 

11 



JATA '/in / 

-SB0XR2 

12 



U=', 

SB0XR2 

13 



I 

SB0XR2 

14 



N 1 0 1> I * c 

SB0XR2 

15 



NTC-G1=3 

SB0XR2 

16 



r.e-3 j 


17 



MC*5 

SB0XR2 

16 



NE=100 

SB0XR2 

19 

c 


MATCH AdCvL numbers TO THE DEFINED DIMENSION IN STORAGE ARRAYS 

SBOXRZ 

20 



ME^Mb.l 

SB0XR2 

21 



MC=2#Mil 

SBOXRZ 

22 



NM=N£+3 

SBi^XRZ 

23 



WRITE (IH,SS) 

SBOXRZ 

24 

c 


READ DATA FGR THE ACTUAL KING 

SBOXRZ 

25 


100 

CALL DATRO(CA) 

SBOXRZ 

Z6 



NEM = 1 

SBOXRZ 

27 



TEST1*0A( 10C1).OA(1002)+OA( 1003)*0A( 1004 1 .DA ( lOOS ) 

SBOXRZ 

28 



CK-OA 1 23 )*0A( 24 1/04(25)46.283185 31 

SBOXRZ 

29 



0I«DAI27) 

SBOXRZ 

30 



L»OI 

SBOXRZ . 

31 



0.^1. 0/01 

SBOXRZ 

32 



DH«0.5*0 

SB0XR2 

33 



WR1TE(IH.60) IDA(I)«I>1,7) 

SBOXRZ 

34 


110 

IF (L) 600,600.120 

SBOXRZ 

35 


120 

IF IM8-L) 600.130.130 

SBOXRZ 

36 


130 

WRITE(1H,65) L.DA424).CK,0A(2S) .0A<23) 

SBOXRZ 

37 



IF (DA<26)) 160.1SO.160 

SBOXRZ 

38 


ISO 

CALL SHAPE 

SBOXRZ 

39 



IFINEM.EQ.2) GC TO IBO 

SBOXRZ 

40 



AC-8.0/AREA ^ 

SBOXRZ 

■ ' 4l 


160 

L1H<H(L<L1 

SBOXRZ 

42 



IF (LIH-HB) 110.170*650 

SBOXRZ 

^3 


170 

LIK2=24MB 

SBOXRZ 

44 



MFLNS » 0A(4<)1 

SBOXXZ 

45 



IFINFLHS.EQ.O} KFLNS ” L 

SBOXRZ 

46 




■ 1 1 HIM 

“%7 



LIH1«2*MB 

SBOXRZ 

48. 



CALL F012(LIHi,LIH2,LP0T,CK,0.A> 

SBOXRZ . - 

49 


180 

CONTINUE 

, SBOXRZ 

50, 



«»0 

SBOXRZ 

51 



Ko>DA(2BI 

SBOXRZ 

52 



CO TO 230 


53 

c 


. * 

SBOXRZ 

54 

c 


PRELIMINARY CALCULATIONS ARE FlNISI^O. 

SBOXRZ 

■ 55 

~c 


THE MEXyTECTIDH IS fcOHE THROUGH TUlTT/rCH MODE." 

SBOXRZ 

36 

c 



SBOXRZ 

57 


'2'00 

“rriOAC26)l 230,210,230" ' ~ " ‘ 

SBOlflfZ 

38 


210 

IF(NEH.E0.21 GG TO 230 

SBOXRZ 

59 


tALL DATRO(OA) ” - • ■ ' feO‘ 


34 



ORIGINAL PAG^ IS 
OF POOR quality 


230 K.-K-J. • 

^ ^ : — = 

IFITESTl.LT^1^01t_CO TO; 250 

«fiITi'UH«55}' ■ ■ 

250 H_RIT£iIWjl5 5_ _ 

HBifiTrH.iibi Ti5Mn,i«i3,i9) ■ ' “T" ~ ■ ” 

15 F0RHAT(1H0*15X«8HM0D£ NO. *13) 

16 Fffft">tAT(lH*»30X*7A10) 

250*280*290 

260 IF (NEK.eO.2) GO TO 240 

Tprint=oa(icoii "■ — - 

call OREt)(Sf OX,SFOy»SFDH,KSFOiSFHX,SFMY*SFMH*KSFH*DA,T 

^ 


SB0XR2 

"5ffOXftZ“ 

seoxRz 

SB0Xk2 

SB0XR2 

■SiraXR2“ 

SB0XR2 


6i 

T5T 

63 

"66 

65 

56 

67 


JITJXRZ 

— ST 

S60XR2 

69 

3F0XR2 

■■ TO 

SBDXR2 

71 

SFOxR?"' 

72 

. S60XR2 

73 

bBQXRZ 

76 

SB0XR2 

75 

SBB)a(2 

76 

S60XR2 

77 

"SBDXR2 

78 

SBQXR2 

■ 79 

5B0XR2 

80 

SBQXR2 

81 

SBQXR2 

82 

SB0XR2 

33 

SB0XR2 

84 

SB0XR2 

85 

SbOXRZ 

■ ' 86 

SB0XR2 

87 

SB0XR2 

88 

S80XR2 

89 

SB0XR2 

90 

SBQXR2 

91 

■3BT3XR7" 

T2 

- SB0XR2 

93 

SBDXR2 

94 

SB0XR2 

95 

SB0XR2 

■ 96 

SB0XR2 

97 

“5BBXR2 — 

■' 98 

SB0XR2 

99 

SB0XR2 

100 

SB0XR2 

101 

'SB0XR2 

102 

SB0XR2 

103 

SU0XR2 * 

164 

SB0XR2 

105 

SB0XR2 

106 

SB0XR2 

107 

S80XR2"- 

108 

SB0XR2 

109 

SB0XR2 

~no 

111 

SB0XR2 

112 

SB0XR2 

113 

SB0XR2 • 

114 

SB0XR2 

115 

SB0XR2 

116 

SB0XR2 

117 

SB0XR2 

118 

G SB0XR2 

119 

SB0XR2 

120 

SB0XR2 

12l 

SB0XR2 

122 

SB0XR2 

123 

SB0XR2 

124 

' SB0XR2 

125 

SB0XR2 

126 

SB0XR2 

127 

SB0XR2 

128 

SB0XR2 

129 

SB0XR2 

130 


”*X X , YY , I W. L , HLC , H; HD", KTBTKmNH ,ii SFS . IPR I HTT' 


Gt>l')«DA't5l) 

290 IFiKJ 310*3 10*3 00 
300 IF ' (rt-ribT*2Cu*jloJ3iO" 
310 CONTINUt 


!F ( Ii- I V( JA ( 50 ) ) .cl.2. a.^D.NEk.EC-1) GO TO 490 

lyO ^"G 

1( I J I 1 1 ) ) “1 1 V, » 32 0 * 4 40 

3 20 

IF?iNT = ji( IC C7T ~ " ~ 

. 4 ALlAA.tr, aToG, YTDG,XTE,SFDX,SFOY,SFOH*KSFD*S*D*CK 
s , 1 -h,l, 4 L, 41 *nc.,nST,FB*HD,HE ."MH *KSFS*TPRINT) 


I r 1 ' (1)1 5'»0, 3oO, 330 

>30 = 1) 

■' CG£‘‘ kilTIO.'^ 

^ . J - 1 , L . ^ t 

■ 5C •( 3 , - ,.) A3 ( J 1 ) .'j/3 . 1415920 5 


>- = ■»,!(*) ^ _ _ 

I 2 A ^2 ( i,, ( ■ *■ - ---- 

!!• - J r-T A- • ( It 1-3 ) 

-.Au- v!.,triAA,tt,Al3G,t7 OG ,xTt,XL£*A, AR * OA H*T *EDG* S* CK*D * YE 

- , lP.il, r, ^LT, '"L— ,AnA 

»-j-;.-*,»_,FL,oL-,,*(l,Fi£F*NST,IK*(i3,l1C*rtE*hD) 


410 


- 20 
22 

24 

25 


I 7 1- '• F I forces FOP LACH MODE 

I* — -■li*l*LJ j ^ ||— L 

•■*17 

••;• : 1 . ( I - V ) (iLtij,r = i,?) ■■ ■ '■ 

.'11- Ln-,'i.',),Ci^,G4(25),OA(23) 

HRITL 1 IW,2C) 


JFlMtH.EO.lJ )ikirLnW.22j 
1FU1LH.E0.2) MRITtlIW,_24) . _ 

hfiITrn»i,25) 

FORMAT 1/lHOlUX* IBHobNERALI ZED FORCES) 

FORMAT I lFl*.28X,22r3 (NO THlCKNEs’S EFFE’CT))’ 
f CRMaTI l,t»,2bX,24H (WITH THICKNESS EFFECT!) 


IPART, IOX,10MA6S. VAL06*6X,11HPFIASE ANGLE) 


00 4/0 H2.1,n 


LALt F ORCU XX,YY.S*SFDX*SFDY*SFOH*KSFO*KSFS*XLE*XTE«YHAX 

f ' ■ ‘ ,AMA,JHAK*HLT,L»NEW,N2,MB,N0*N8*NM*CK*T) 


ilAAC*I(l.l) 

52=AC4T( 2,1 ) 

33= 50«T(S1*4a+S 2^F2) 
i 4*5 7.29570* AT AN2( 52. Si) 


35 



n n o 


<i50 

30 

<i70 

480 

490 

HRIT£IJH,30) f<l,H2,Sl.S2«S3»S4 

FCPnA-.T (lh02le*lP3£l'9.5,0PlFl6.4) _ . • _ 

LUMl-NOt — - . . - . . . 

mkiu nw,55i 

SB0X<^2 

SB0XR2 

SB0XR2 

SB0XR2 

131 

132 

133 

134 

135 

136 

CDNMNUE 

1M0A(26).G1.C.) GO TO 510 

SB0XR2 

S8DXR2 


OO 500 I>13«22 

SB0XR2 

137 

&00 

OA(I)“7 

SB0XR2 . 

136 

510 

If lIFIXIDA<5Cn-l.E0.0) GO TO 100 

SB0XR2 

139 


lFtNEh.tQ.2) GL TO-lOO 

S80XR2 

140 


NtH»i 

SB0XR2 

141 


tF(uA(26>] 18Ci5aOfl80 

SB0XR2 

142 

580 

IPR!NT-DAiiC04) 



call MKLD(DAtT«NHyNBtKSFH,SFHX«SFHY«SFHH« lUflPRlNTl' 

SB0XR2 

144 


CC TO 150 

S60XR2 

145 



SB0XR2 

146 


iKKilR EXITS 

SBOXR? ■ 

147 



SB0XR2 

148. 

590 

I PR-30 

S80XR2 



GC TO 610 

SB0XR2 

150 

600 

I PR-2 7 ■ • - " - ~ 

Sb5XR2' 

151 

610- 

WRIT{.(Ih,J3» IPk 

SBQXR2 

152 

35 

FUK.rtAT ( lh010X,eHBAD DATA14) 

SB0XR2 

153 


GO TO 700 

SB0XR2 . 

154 

650 

KRITt (!H,4C1 

S80XR2 

T5T 

40 

PClK«AT(lh010X,42HLATERAL LIMIT ON NUMBER OF BOXES EXCEEDED.) 

S80XR2 

156 

700 

5T0? 

SB0XR2 

15? 

55 

FOHMATdHU 

SB0XR2 

158 

60 

FCRMATtlH0,10X,7A10) 

SB0XR2 

159 

65 

FCRMAT(lHCtl0XiI2i23H BOXES ALONG ROOT CHORD* 

SB0XR2 

160 

I 

15X,Z2riR00T CHORD LENGTH -*F8.2>3H FT/ 

sbUxRz 

161 

i 

lHO,iOX»19HKtOUCED FREQUENCY -*E6.3> 

S60XR2 

162 

3 

' 15X,Z2HFR'eE STREAM VELOCITY -tF8,2*7H FT/SEC/ 

SB0XR2 

163 

4 

lriC,10X,llHFREOU£NCY - , 1 P Ell . 3, 1 OH CYCLE/SEC) 

SB0XR2 

164 


INC 

SB0XR2 

165 


36 



ORIGINAL BAGD IS 
OF POOR aUiilTY 




.SVBRO^qHK' OATKD(OATA) 

OATRO 

2 



Kiftt-REAfJ 5UBROU11NE -fiAtROf OATAt 1 1 )" 

DATRO 

— -j- 



0Ui£ltS2fiM lmS>Jtl2)«l>ATA(l)t00KBU(I0) 

OATRO 

4 



BATJT -A1^1^&HAU>HA/«0TEST/1H /,ETeST^lH-/«SIEST/lH*/ 

C^ffTBiir 

5 



DATA mu / 

OATRO 

t> 



l*«5 

OATRO 

~ 7 



lW-6 

OATRO 

8 


1 

■rSaS fi%«2) EHlH^ALPtlNOflORBUilltl'l*?) 

OATRO 

9 


2 FeRMA'HAl«A5tIfa«6A10f AS) 

OATRO 

10 



IFUHO.EQ.O) CO TO 20 

oaTrd 

11 



IF tEHIH.NE.SIESTt GO TO 10^ 

OATRO . 

12 

~ 

" " 

NEU SING IF COLUMN 1 CONTAINS A PLUS SIGN 

OATRO' ■ 

13 

c 


1NITIALI2ATJQN Of DATA ARRAY 

OATRD 

14 



00 l6i I°23«1CC3 

(TiTRO 

15 


101 

DATA! Il-O.O 

OATRD 

16 

— 

~~ " 

DO 102 I>1.22 

OATRD • 

r? 


102 

0ATAU)>2 

OATRD 

18 

, - - 

■“ 

OG 103 I=>=10A, 7C0,A 

■ • OATrO- 

. 10 


103 

OATAlD'l.O 

OATRD 

20 


TOT 

thwttnue “ ■ 

- -jj-JfTRD 

~ZT 



IF (ALP.fcO. ATEST) uO TO >5 

DATRO 

22 



IF tALP.Nt.CTESTJ 1,3 TO 3 

OATRO 

23 

c 


numeric carc 

OATRD 

24 



CO 3‘r=i»F 

OATRD 

25 



uCKBUI I ) =of eU 1 I J 

OATRO 

26 



LCM IhiJt 

DATRO' “ 

27 




OATRO - 

28 



L,: 3 1 = .,; 

OATRD 

29 



1 F 1 

DATRO 

30 

c 


FJk 

DATRO 

31 


- 

i f < S I sjt"* 

OATRD 

32 



f T i 1 i n;, j = , , > L U ) 

DATRD 

33 



= ;.■>. 

OATRO 

34 



T ; 11 

DATRD 

35 




OATRO 

36 



r : : ' i - 1 , ? 

DATRD 

37 



_aTai i ■>*„) = , - 1 li i 

. OATRD 

36 



- 1 

DATRD 

39 


I 1 

IF I-.'Ih.Ni.eIijT) j, T.‘ 1 

OATRO 

40 

c 


I' i _:-<TAlNS A fl.NljS >ICN 

DATRD 

41 



- z "t Jt '' 

OATRD 

42 

• 


t‘.’ 1- -ii-- 

OATRD 

43 



p-»ir 

OATRD 

44 



.iE- -UT 

DATRO 

45 



S T ^ ' 

DATRO 

46 



1 A V ^ 

DATRD 

47 



C ^ 

DATRD 

48 



^ r 1 T t 

DATRD 

49 



. r if-. , 1 N j < 1 3Rr V, 1 1 1 * I = l > 7 ) 

OATRO 

50 




■ DaTRd 

“5T 



s : ^ . 

DATRO 

52 


4 

F . ► --T ( -E i:.E 1 

OATRO 

53 


4 ; ; 

r T 1 . -•• -M AIa „■< T-li Ei..b. JOi TERMINATED -) 

DATRD 

54 


-f j 

► I'lAoi - ►Ai !<■»» ?AiO) 

DATRD 

55 


-* 

f„'-iTt///lh tli'/tt/"' ""i-: r^ATACA^D/) 4 

DATRO 

56 



t 

“■ “ DATRO 

5T 


37 



KST 

IBtNH 

.LtNEW 


1 . 

I l ... i 
l'clIt = lP‘'i.'i/i.C!i. 

Jtr I''=^1 PkI N- 1 L ‘..O*-'*'- ! 

Jk m «!P«lN/iCC 

: 

Ik lTiilP^!ki/lC ' 

iMNtw.ea.i) iRnt = i_PxiN-iJ«^ipiTc ^ , . 

If JNlH CC TO ^ j^ 

lLDb = C „ 

ns=D4(.;'i ) 

IF <Nb» abo.oSCtioo 

100 "if {NS-NEJoi*n lo^riOsVtfSo 

10-: NiP»f^*i 
NS=N$P' 

IF 1JA (24)I O55»b53.110 _ 

no 00 115 I=l,p.iF 

X£OGtH=SA(a»I*27>/OA(2^) 

115 Y60GU l=0A<2*I + 26)/DA(24) 

XED GI1) = 0.0 

HSf-0AIA5) 

IF (NST» 871» 125.> 130 

'125 XTDG(1)>1.0 

NST»i 


GO TO 141 

1 130 IF <WST-MTDGi*lJ 135-tl35»871 


135 NSIP-HST 

DO 140 1*1»NSTP 


YT0CIII*0A(24I*43)/DA(24> 
140 XTOetll>OA(2*I«44)/OA(24) 


141 «ST»«H5T+1 
M5T*MSTf 


XT0CFHSTP1-XEDGIH5PF 
YTOCtWSTPI»VECC<HSP> 

VTa&(iJ«o.o 

IF <HST.£Q.MTCGI> GO TO 150 


JL^NST*-! 

00 142 t*A.*HTOGI 


XT06111*XTOG(I-l) 

I 142 rT06tII-YT0£Il-l> 


150 IFfASSOrtOGIHSTP >-XEOG{MSP MYTD6INSTP I^YEOCINSPl I .CT.ERRRI 

1 ■' $ GO TO 871 


FQR'ifING WITH THICKNESS 


SHAPE 

A 

A 

A 

A 

A 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

sha'pe" 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

SJO^E 

■ shape" 

SHAPE 

SHAPE 

SHA^ 

SHAPE' 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 


SHAPE 


KAPE 
SHAPE 


SHAPE 

SHAPE 


SHA 
SHAPE 


SHAPE 

SHAPE 

SHAPE 

SHAPE 


SHAPE 


H 

SHAPE 


SHAPE 

SHAPE 


SHAPE 


HA 
SHAPE 



DO 215 1=1, L 

' IF jMLC<l,l).fct.0J GO TO 215 
J*hLC(l»U 
K=MLCC2,U 

CALL SURP 2(XXU) ,YY, J,<, 1,AHAC1,1 t ,DM«DN,SF)iXt SFMY^SFWH .KSFH, 1 1 
215 CONTINUE " 

XMJKiIt.EG.ul GO lu 235 ; 

PRINT CUT NflCh DISTRIBUTION 
KRITE(IWffeO) 

DO 230'1-1,L 

IF XNLCIJ,IKEC.qj GO TO 230 


SHAPE 


yiAPJ_ 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

“SHAPE" 

SHAPE 

SHAPE 

SHAPE 


38 



jn o 


ORIGINAL EAGi; IS 
OP POOR QUALITV . 


4L»HLC<2,J»-MLCa,U*l ' 

IFUL.Ea.Oi GO TO 230 

weiTEIlMtTOI I 

JLP^L/fe _ ■ - _ 

If fjL-6«JLP-f*e.O) JLP=JLP*1 

Kl»i1LCtIil 1 -1 

JL-m.C*2,I) 

DO Z25 J«1,JLP 

K^nrm 

225 HRITf « I«,805 t I J 1 , AHA I Ji , I ) 1^J1«K • JL, JLP ) 

zao'cGNfihuE 

-Zas'CdNTiNUE ■ ' 

&U lu 7bo 


SHAPE 
SHAPE 
SHAPE 
SHAPE 
SHAPE" 
SHA>£ 
SHAPE 
SHAPE 
“SHAPE 
SHAPE 
SHAPE 
SHAPE 
"SKA PE 
SHAPE 


YWAX=»rh DC«HSP) 

TFIVl.LT.OT?n"“{ 

AREA'0.0 

XX«l»»Dri 

YXCU-OH 


~'Td 'dbO 


SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 


YYU)-YY«I-1» + D 


•CAL£uTATE“PJ?TriaL'FUOr-¥lTlG-ART5- AHB~ firTNEM,I ) 

COnPUTE AND ICENTlfY 80X OISTRIBUTIQH QN THE MING PLAHFORH 

■CaCl PLNEH(xEOG,y£DG«xTOG,Vn)G>yy.VViiXLE,XlE,AR;A'KkA, 

1 ^MJ^LC >MLW, HEH» HS »NST<1R1 TE> 

If IMLfLl.GT.Pe) GO TO 825 


SHAPE 
• SHAPE ■ 
SHAPE 
■'SHAPE 
SHAPE 
“SHAPE 
SHAPE 
■ SHAEF 
SHAPE 

“SHAPr 

SHAPE 


00 570 l>ltL 
J«AK° HAXO( JHAKyMLCI 2< I ) » 
570 JHAX-HAXOt jHAXtHL (I)> 
JHX^JHAK 


SHAPE 
“SHAPE ■ 
SHAPE 
SHiFE" 
SHAPE 


00 bZO J«1«JHAX 

EDGt J> = 0.0 

BO 620 K»l»2 
HLT1K»J1=0 
11 = 1 

1 

Ki=a 

OC o55 J»ljJMX 

_IF tHLC tZ.Il l .Cfc. J ) G0XU_632_ _ _ 

Il»[fVl 

GO TO 630 _ _ ___ _ 

HZl=ll’ 

flHO CRBER Of trailing EDGE bOX AT J-Th CHORDWI S£ ROW 
IFtIZ.GT.L) GC TO o38 
If 1HLC12,12) .CT.0)_K1=1 
I'F t<i.E0'.0"l“GC fo' blV 

IF MLCl^,!.^) ._EO.O.ANO.NZI.N£.OJ[^ GO TO 65b 
KLS = I A&StHLHtlZn 

<HL>HI2>.NE.0) go TO 63b 

I2=I2^i 

K2 = i _ _ _ 

GC to b35 

IFtHLS-JJ fc 55.fc 5b,o38 

IF t'ki.LT’.O) GC to 053 

I F tHLWUZl ) t 5Z.65Z i65b 

I2=I2*hZ 

iFtI2.L€.L» GC TO 635 
N2Z=L 

IZ • 

CC"'tO 65z' ‘ 

n2=MLmi2)yp.LS 


SHAP E 

SHAPE 

SHAPE 

SHAPE' 


SHAPE 

109 

■“SHAPE 

■ 116 

SHAPE 

111 

SHAPE 

• 112 

' SHAPE ■ 

113 

SHAPE 

llA 

SHAPE 

115 

SHAPE 

116 

'shape 

' 117 

SHAPE 

118 

SHAPE ■ 

" 119 

shape 

120 

SHAPE 

121 

SHAPE 

122 

• SHAPE 

123 

shape 

124 

SHAPE 

125 

SHAPE 

126 

SHAPE 

127 

SHAPE 

128 

SHAPE 

129 

SHAPE 

130 

SHAPE" 

131 

SHAPE 

132 

SHAPE 

133 

SHAPE 

134 

SHAPE 

135 

SHAPE 

136- 


39 



GO TO oA2 
6^6 NZ2«^2-1 
’GO'TO'biTi 
b^ia N22=I2, 
b^Z CCitTi.noc. 

HLTU,J)»NZ1 


SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

137 

138 

■ 139 
160- 

■ 161 
162 

■ «LT12»J)=NZ£ 


SHAPE 

■16 J 

655 CQNTINOE 


SHAPE 

166 

C 


SHAPE 

165 

C FIND ORDER Of RING bUXES AT EACH I-TH SPANRISE COLUMN 


SHAPE 

166 

DC b70 i=i.L 


SHAPE 

167 

IF t.RLHm) b6'%,662»b62 


SHAPE 

168 

662 HLCC1,I)^MLH(1)«I 


SHAPE 

169 

IF IMLClZf U.EC.OI Hi.C(l.n«0 


SHAPE 

150 

GO TO 670 


SHAPE 

.151 

66^ NLC(2,n-IABSt«LRm)-l 


SHAPE 

152 

MLC{I,II-=l 


SHAPE 

153 

670 CONTINUE 


SHAPE 

156 

C 


^A^ 

T5F 

IFURITE.EQ.O) GO TO 678 


SHAPE 

156 

MRIT£(IW,b5) 


SHAPE 

157 

IF(NEH.£Q.U hRITEtIR,77) 


SHAPE 

158 

IF<NEH.E0.2) RR1TE(IH»78) 


SHAPE 

159 

WRITclIW,e2) 


SHAPE 

160 

00 672 


SHAPE 

■uai 

672- HR1TEIIH«85) J«HLT(lf Jt tHLriZ* J > 


SHAPE 

162 

MRITE<1H*83) 


-SHAPE 

163 

00 676 J«1»L 


SHAPE 

166 

C 676 URITEIIH.85) J«FUJCU.*->>«m>CI2« Jl 


SHAPE 

165 

676'HRITE(1R«B5> Jl ,HLC{2«J1*HLM< J1 ,NL( J) 


SHAPE 

166 



SHAPE 

167 

e ^ 


SHAPE 

168 

IF(XE0O(NSP».CT.i.l GJ TO 365 


. SHAPE 

169 




SHAPE 

170 

DY = Y£DGlNSP)-YiL!ClNirt) 


SHAPE 

171 

IF lAdStOYI .LE .LRR^) IEJG=I 


SHAPE 

172 

IFIOY.LT.O.I CC TJ 627 


SHAPE 

173 

C CCHPuft VALUES FGk leading EOgE CCKRECTIDN 


SHAPE 

176 

YMAX2=YrtAX">YNA> 


SHAPE 

175 

DY=0 


SHAPE 

.176 

DC 72C J = l,Jrl> 


SHAPE 

177 

IF (lEOG) 710,715,710 


SHAPE 

178 

C OY If) FOLlCRING tXpKzSilON IS AkdlTKA»Y 


SHAPE 

179 

710 ECGI J)=SORT( (YrAX2-rY( J )*YY( J) ) / 1 D Y« ( 2 , A YUA X-0 Y ) I ) 


SHAPE 

180 

IF <OjG< J.I.CT.l.u) ;0 TJ 71; 


SHAPE 

181 

GO TO 720 


SHAPE 

182 

7 15 L DGI J )- 1 .G 


SHAPE 

183 

720 CCNTIhbt 


SHAPE 

186 

C 


SHAPE 

185 

750 CCHTINbt 


SHAPE 

186 

C 


SHAPE 

187 

7W CDfiriNUt ■ ■ ‘ • ' 


SHAPE 

188 

c 


SHAPE 

189 

RETURN 


SHAPE 

190 

C 


SHAPE 

191 

800 ILFRji< = 800 


SHAPE 

192 

T( 1,1)=X£0GIK) 


SHAPE 

19 3 

H2,U=YcDGIK) ■ ■ ------ 


SHAPE 

196- 

!ER = : 


SHAPE 

■ 195 

GTi T 3 i5 ^ ^ 


SHAPE 

196 

305 l£PRQR=a05' 


SHAPE 

'197 

GC TO 860 ■ ■ " 


shape 

198 

810 IERK0ft=8i0 


SHAPE 

199 

T ( i T 1 ) - UA ( ^ / i 


■ SffXPt“~ 

~ Tinr 

T(2„1)*NSP 


SHAPE 

201 

I£T^= '2' - ■ • - ' -- ■ - - 


SHAPE' 

20^ 

GO TO -860 


SHAPE 

203 

815 le'RKQR^grr ' ■ 


SHAPE 


GC TO 860 


SHAPE 

205 

820 lEk-^z ;rrz ■ — 


SHAPE 


T(1,1)=X£0G(NSP) 


SHAPE 

207 


ko 


OKIGINAL EAGi) Ig 
OP POOR aUAMTy 


T! 2 , 1 )-YEDG(NSPT 

T! 3 *i)*NSP 


SH*PE~ 

SHAPE 

~zor 

209 

niT^pX ■ 

GO TO 840 


“ " SHAPE 

SHAPE 

zie 

211 

“Tf 2 b‘ rERftUR«tt 25 



■ "ZIZ 

TI 1 -, 11 =nEw 


SHAPE 

213 

TT 2 ,Iii*“ - - . . . 


SHAPE 


Tt 3 tl)«HL!K) 


SHAPE 

215 

t< 4 fl)«ML!K-l) 


^ SHIITE 

216 

IER =»4 


SHAPE 

217 

■ GO rOTVO 


SHAPE 

Elf" 

627 lERR 0 fi >627 


StttfE 

219 

T( 1 » 1 )«>TE 0 C!NSP) 


SHAPE 

Z 20 

T( 2 . 1 )-YEDC 1 HSI(» 


SHAPE- 

221 

T]I 3»1 )>NSP 



~ZET 

T< 4 « 1 )-HSH 


SHAPE 

223 




~»r 

60 TO 640 



223 

630 I£ft{tOR >736 




T(lfll>I 




TI 2 « 11 -J 


- V 


T( 3 * 1 T»X 


-5 


T< 4 ,ll-Y 



- 

TI 5 fll-EH 




T( 6 . 11 >£MY 


. 


T( 7 « 1 I* 0 EL 



• 2 ^ 

IER -7 



239 

640 MRlTE!IH* 20 t lERROR* CT 4 I « 1 », I« 1 *IER) 



23 ® 

STOP 




850 IPR -29 




GO TO 890 


shape 

239 

655 1 PR «24 


shape 

239 

GO TO 690 


»(APE 

296 

360 IPR >30 


SHAPE 

291 

GO TO 690 


SHAPE 

292 

865 K> MIN 0 (K»NS) 


SHAPE 

243 

IPR« 2 4 K +29 


SHAPE 

299 

GO TO 890 


SHAPE 

295 

871 IPfc -45 


■ SHAPE 

296 

690 NRITE 


SfflTEE - 

29 ? 

0 nw.ioiiPR 


SHAPE 

298 


10 

zo 


FORMAT! 1 HO, 10 X, 17 HSHAPE — BAD DATA«IS 1 
FORMAT! 1 H 0 » 1 CX> 38 HBAD MOHBER IW SHAPE HEAR STATEMEHT M 0 .»I 5 « 


SHAPE 

SHAPE 


2A9 

250 


1 /IH ,l&X,lP8fcl4.6J 

40 FORHA T!/l HOtSX<3HNO.Tl2 ,43H R ED 1-STRlBUT IQH OF RIHG LEADING EDGES* 
■ ■ 'i NSI.li»4)iJ - iI2»7H, lEDGi3H « *111 
45 FO RMAT(/lMO»SXiZ ^HH I MG TR A ILlWG EDGES ) WST I t11»4H ) 

»■ lOh, ItDG »— ili> 

50 FCR«AT( 4 C 6 X.I 3 ,lP 2 £ll, 4 n 


SHAPE 

SHAPE 


251 

252 


iI2« 


SHAPE 

SHAPE 

SHAPE' 

SHAPE 


253 

254 

255 

256 


60 

65 

70 

75 


FORMAT! ittl.SX . 47 HL 0 CAL 
FORMAT! IHl) 


MACH HUMBER DISTRIBUTION ON PHYSICAL 


FGRMAT! 1 H 0 , 5 X, I 2 , 19 H-TH SPANHISE COLUMN) 
F ORMAT! 1 HC^ 6 X ,I 2 t 1 7 h-TH Ch^ROHlSE ROM) 

77 F 0 RrtAT!' 4 OX, 13 F>H¥S'lCAL“HING// 7 ) 

78 FCRMAT 140 X 1 16 FTRANSF 0 RME 0 MING///) 


and LAST! TRAILING) 


80 FCRMAT!1H , 5X , 6 ! 2X, I 3 »1PE13. 5 > ) 

82 FORMA T!15 X,6ehO RDER OF F IRST 1 LE AD INC) 

IIN CHOROn’ISE RCwV/ 2UX« IHJ, 3X, 12HMLTIN 
33 FORMAT! // /15X ,fe3H0R0ER OF FIRSTIRUCT) AND LASTITIP) KING BOX 
lANnI5t~CbluM)r/72o‘x,lrilT3'x,12HMLC(hEW,l, I ) , 3X , 12HMLC! NEN 1 2 ♦ I ) 
2 |4X<llH(MLMStHti)v5X>10HML(NEht]) )/ ) ■ 


35 FCfi«AT!l6X,I5,5X,I5i J!iOX,I5)) 

86 FORMAT! 20X.56FCAUSSIAN INTEGRATION 


POINTS IN CHQKOMlSt ROM - 


IM,K,31//I 

37 FORMAT! ///2Ca ,52HGAU3S_1 AN 

'ifKI/V) 

B8 F0RmAT! ///20X,60h5PANMISE __ 
IT GX!2iAiJ)//I 

89 FDRMAT!1H i 5X_j 2 1 2 X •. I J , 1 P2E 16- 6 ) > 
STOP 
End 


Integration POINTS IN SEmI-span - 
CUURJINATE (GYP) AND MACH NUMBER 


HING/) 

SHAPE 
SHAPE 
'SHAPE ■ 
SHAPE 
SHAPE 
SHAPE 

257 

256 

259 

260 
261 
262 


SHAPE 

■'263 

C BOX 

SHAPE 

264 

,J)/) 

SHAPE 

265 

. IN SP 

shape 

266 


SHAPE 

267 


SHAPE 

268 


SHAP^E' 

269 

GX!NE 

SHAPE 

270 


SHAPE 

271 

GY!NEM 

SHAPE 

272 


SHAPE 

273 

CMP ) A 

SHAPE 

274 


SHAPE 

275 


SHAPE 

276 


SHAPE 

277 


SHAPE 

278 





iilM 


OIHENSIQN ,XfcB6<lii»»£OGrSr*X?B6< 


e*TA Zl/iHf«/»Z£/2HG«/ 
i»M>.5*0 




NSP «N5 
. NSTP»I»5T 
00 30 1>1»M8 
MUf2fll«0 

HLtin ) = 0 


XLE C 1 )»X 
XTEm^^XEOGtNSP) 

00 30 j>i*tie 

30 AR(Jtl)«0.0 


210 Kl-0 


YR«YE06U> 

X2-XR 

Y2=YR 


XE 

sun>o.o 

-“ar-ex^^to — 

ASHG<=1,0 


PLIV<f 




PLKFM 

«J<FH ■ 


PLNFM 

"Fcranr 

PLI^H- 
■~»LHfh ■ 
PLNPH 


PLNFH 

T^TiFS 

PLMF H 

PLNFIi 

PLHFM 




PLMFIt 

PLMFtf 

PLNFH 

plsfh 

PLNFH 


rLwn 

PLNFH 

PLNFH" 

PLNFH 

■ptWH 

PLNFH 



ICH0=0 

"IFIN=0 

IE0Z=0 

XL“XR 

YL-YR 


r< * ■ s n IT 5 1 


IF(ASNG) 225«230»230 

"225 IFIXR.lt. XTOG fl). AND. XTDGiiJ-XR.Gt.ERR2) GO TO 530 
IF(K2.EQ.01 CG TO ^40 


GO TO 

_^30 CQ HTIHUE 

IF ( K 1 ♦ 1 .G I . N S I GO" “to” 4 2"0 
IF(Kl.EO.O) GC TO **00 
240~ICHO=0 

IF(XR.GT.XE2) ICHQ=1 


JflN-0 


Y=YYl Jl+Ori 

IFCASNG) 252.256.256 

252 lF(ABSfXTDG(i]-XTOGCNSTPI).GT.ERR2.0R.I.NE.L} GO TO 254 
IF(Y-DH.GT.YTGG<NSTP 1 1 GO TO 530 ■ ' - 


XT£«J)*XT0G(1) 

GO TO 250 ; ' 

254 SUH~SOH-ARIJ-,U J ' 

GO TO 258 

256 AR<J.1)=0«0 

258 CONTINUE 

IF(YL.L T.Y ) GG TO 2W _ 

IF(A5NC) 265,270,270 -- - -- - 

265 AR1J.I)-0.0 ^ 

GO TO 250- 

270 CO NTT*<UE _ 

ARtJ.lJ-i.O’" 

SAH=SAH *AR< J,n 

CO TO 250 


PLNFH 
■ PLNFHT" 
PLNFH 
■PCWTT 
PLNFH 


PLNFH 

FCHFIT" 

PLNFH 

PTNFH' 

PLNFH 


runrn 
PLNFH 
PTHFIT' 
PLNFH 
“ mSFfl' 
PLNFH 



PLNFH 

54 

PCHFiT 
- PLNFH 

5S. 

56 

“PORT* 

57 

PLNFH 

58 

PLNFH • ■ 

50 

PLNFH 

60 

PLNFH 

61 

PLNFH 

62 

PLNFH 

63 

PLNFH 

64 

PLNFH 

65 

PLNFH 

66 

PLNFH 

67 

PLNFH 

68 

PLNFH 

69 

PLNFH 

70 



of fee® IS 


260 

IFtV-Yt-*ERR2.6£.0» CO TO 290 

PLNFH 

71 


AR4 Jvl)>«R( J« I)«ASNC*(rL-T40l/0 

PLHFN 

72 


Xl»Xt. 

PLNFH 

73 


Y1«YL 

PLHFN 

74 


1F( ICHC.NL.l) GO TO 310 

PLHFN 

75 


iFOKl 284,. 85,285 

PLHFN . 

76 

284 

IFiY€2.LE.r-0 J GO TO 320 

PLHFN 

77 


GO TO 286 

PLHfH 

78 

285 

[FtYE2.GE.Y} GO TO 320 

PLHFH 

79 

286 

COHTINUE 

PLHFN 

80 


IF(XE1,CT.Xl) go to 400 

PLHFH 

81 


IE0Z=2 

PLNFH 

82 


GO TO 400 

PLHFH 

83 

290 

X1-X2 

PLHFH 

84 


Yl«Y2 

PLNFH 

85 

295 

IEDZ-0 

PLNFH 

86 


1F( IChC.Hl. 1) GO TO 310 

PLNFH 

87 


YR»Yt2 

PLHFH 

88 


GO TO 315 

PLNFH 

89 

310 

YR«Ytl4TNC*«XR-X£l) 

PLNFH 

90 

315 

1F<JX) 317,316,316 

PLHFH 

91 

316 

IF(YR.GE.Y) GC TC 320 

PLHFH 

92 


GO TO 318 

PLNFH 

93 

317 

IFIYR.LE.Y-O) GO TO 320 

PLNFH 

94 

318 

CONTINUE 

PLHFH 

95 


[FUCHC..EO.U GO TO 325 

"PLWH 

96 


IF1N«1 

PLHFH 

97 


X2-XR 

PLNFH 

98 


Y2»YR 

PLNFH 

99 


GO TO 330 

PLHFH 

100 

320 

JFIN=1 

PLNFH 

101 


IF(JK) 322,321,321 - 

plhFh 

Id 2 

321 

Y2=Y 

PLNFH 

103 


GO TO 323 

PLHFH 

104 

322 

Y2=Y-0 

PLHFH 

105 

323 

X2=XEl*lY2-YEi)/TNO 

PLHFH 

106 


IFeA8SIX2-XRI.GT.ERRR) GO TO 330 

PLHFH 

107 


If IN-1 

^h¥h 

l08 


JFIN-0 

PLNFH 

109 


GO TO 330 

PLHFH 

110 

325 

lEDZ-1 

PLHFH 

111 


X2-XE2 

PLNFH 

112 


Y2-YE2 

PLHFH 

U3 

330 

AR<J,I)»ARe J,I )♦0.5*ASNG4^2.♦XR-X2-Xl)•^y2-Yl»/02 

PLNFH 

114 

c 


PLHFH 

115 

c 


PLNFH 

116 


I f t A:? ( T NG Ti'iG=i*t — cO 

PLNFH 

117 

1 331 

CONTINUE 

PLNFH 

118 


IFeJ4) 332-,J31,3 j1 

PLNFH 

119 

331 

iFfJ.LT.JLT) CC lu 334 

PLNFH 

120 


xLT-Atl*lYY(JLn-Ycn/TNG 

PLNFH 

IZl 


IF<XLT.GI.Xt2.*tvj.ASNG.OT.O.O) oO TO 334 

PLNFH 

122 


If (YtOGe i7.GT.C.0.Ain).yY( JL7 1 .Ll .YL. AKIr.ASNG.GT.O.O) XLT-0.0 

PLHFn 

123 


If (YY( JlT).uE .YTIjGCrtirpn xn uGINSTP ) 

PLNFH 

124 


If (ASNG.CT .C.C) XL=IJLT)«*>.T 

PLNFH 

125 


IF (ASNG.lt. G.C) aIlIJlTJ-XLT 

PLNFH 

126 


oC TJ 3 33 ■ • ■ 

PLNFH 

127 

332 


PLHFH 

128 


XLT-Xtl*nY ( J,.T )-Yt n/Tr«0 

PLNFH 

129 


IFeVY<jLT).GL.Yli,o(NjI^H ALl = AlGGlNSTh>) 

PLNFH 

130 


XTEtJLT)=XGT 

PLNFH 

131 

333 

JLT-JL7 ♦JK 

PLNFH 

132 


GO TO 133i 

PLNFH 

133 

334 

CCNTiNUt 

PLNFH 

134 

C 


PLNFH 

135 


If UiLiZ .lG.I ) CO ru 4)0 

PLNFH 

136 


iAn = 56N,ASI j, 1 ) 

PLNFH . 

137 


IF(Jk) 344,34£,342 

PLNFH 

138 

-u?- 

rCNT'lNUE 

PLNFH 

139 



PLHFH 

140 




■ ' -GlTTj'35'6 • ' 

■PLHFH 

141 


3A<i IMJUN.EC.i) CC Tj A ?t> 

PLNFN 

142 


3b0 IFtASv^J 2 s 5J i260,3A0 ■ 

PLNFM 

143 


3^5 IF 1 If iN.fc.l ) GO lu Hli ■ ■ 

PLNFN 

144 


CC TO o30" ■ . ■ . .. 

PLNFM 

145 


360 IF UFlN.EO.l > CO T J A66 

PLNFN 

146 


GD IJ 600 ■ 

■FCHFH ■ 

147 

c 


PLNFN 

' 148 

c 

LiA)!NC tt-Gt TftAlulNG EDGE CONPUTATIOK 

fiLMFN 

149 

c 


PLNFN 

. 150 


AOO ifiasmCT aac»a’C5,a05 

PLNFM 

151 


A06 IFiJK) ABOtAlCtAIO 

PtNFN 

152 


"ATTr K1-K1*JK — ■ 


153 


IFCKl.GE.NbP) GO TO A<!0 

PLNFN 

154 


jrer^Cinmni ' - ‘ 

PtNFIT" 

isy 


YEl=YtOGIKl) 

PLNFN 

156 


K-K1*JK 

ULTIFir 

' 157 


XE2=XfcOG(iU 

PLNFN 

i58 


YF?=ytuc(Tn 

“PClffN 

157“ 


GO TO A50 

PLNFN . 

160 


A 20 XF'fjK.LT.O) ■GC"W"jA(r' 

“PCIBFfl — 

I6T' 


JK«-l 

PLNFN 

162 


JLt«J 

"~FLWH 

163“ 


KLCT^*iJ=0 

PLNFN 

164 


K2=K2-JK 

FUWH 

I65T 


430 K2-K2+J>< 

PLNFN 

166 


XEl=XTDG(K2l 

PLHFH — 

167“ 


YE1-YTDGCK2) 

PLNFN 

168 


K«K2*JK 




XEZ-XTDGIK) 

PUIFN • 

170 



PLNFN 

171 


GO TO 450 

PLNFN 

172 


440 K2*K2+i 

PLNFN 

173 


XE1-XTDGIK23 

“POSH 

176" 


Yfcl=YTDGtK2I 

PLNFN 

175 


XE2«XT0G(K2*11 

PLNFM 

176 


YE2=YT0G1K2+1» 

PLNFM 

177 


450 6-XE2-XE1 

PLNFM 

170 


F-YE2-YEI 

PLNFN 

179 


TMG-F/IG*l.t-20» 

PLNFN 

180 


IFtG.LT.O.} CC TO 650 

pLhfh 

ISl 


AREA»AREA*ASMG«G#IYE2*Y£U 

PLNFM 

182 

c 


n.NFH 

183 


ICHQ=0 

PLNFN 

184 


IFtXR.GT.XE2) 1CK0»1 

PUffN 

109 


IFtlEOZ.EO.L) GO TO 290 

PLNFN 

166 


1FI1E0Z.EO.2) CO TO 295 

PLt^N 

187 


GO TO 250 

PLNFN 

188 

c 


PLNFN 


c 

number OF BOXES IN SPANWISE COLUMN 

PLI^H 

190 

c 

HAKE TYPE AND ORDER OF HAKE BOX 

PLNFN 

191 

c 



192 


465 IFfASNG) 478»467t467 

PLNFN 

193 


467 IFtJK) 476»46S»469 

PLNFN 

194 


469 JJ-0 

PLNFN 

195 


IF tARt J.II.LT.ERR2) JJ-1 . . 

PLNFN 

196 



PLNFN 

197 

, C 


PLNFN 



JJ«NL(I) 

PLHFH 

199 

1 

DO 471 

PLNFN 

200 

J 

KK“K 

PLNFN 

201 

1 

IF tARCK»I).LT-.AHAKl GO TO 472 

PLNFN 

202 


471 CCNTINUE 

PLNFN 

203 

' 

HLCt2.H=Ki< 

PLNFN 

204 


GO TO 473 

PLNFM 

205 


472 HLC(2»i‘)»Kk-l 

PLNFN 

206 


473 CONTINUE 

PLNFN 

207 

C ADJuSTl1t^T iii- flRST dOX AT Lt^AOING EDGE ALONG KK-TK CHOROHISE ROW 

PLNFN 

208 


KK’HLCI 2« n 

PLNFN 

209 

•* 

IF (XXCn.Gc.XLE (KK)) GO TO 474 

PLNFN 

210 



O U W 


ORIGINAL BAGE IS 
,OE POOR quality: 


r>tCt2 il MML CtZ<l l-i. 

“CC ' M2 
CONTINUE 
r 

U 1 u 7 C ii 
476 nLS*-l 

NH I ) »ML( T-1 ) 

K=NL«I) 

G L" T 3 4 d o 
473 NLb=l 

4S C 'CG 47 V J J= i', ■< 

Ir (AKCJJtll) 4d 1 j‘id' 1 , 4 04 

431 IP t A oSTaP IJ J ,'U )-cK 4 2 ) 4b3.4dit4e2 

C _ 

432 if? ( InTSo) Ekx2, JJ,I ,A^IJJ,U 


"433 A<!lJJ,n=C.L 
454 CCNTiNCt 

Ir t-L«n J .Ni .C) «U 

iPCAiNoT 4 o 7 > 4 3'b» 4o6 
C -IC UP ^Ah-e. 

400 Ir ( A A i JJ , I > . Cl .ip<A< ) GO TO 490 

"L"ti>AJJ _ 

G'J T' 4oS ■ 

C "iNr rJCT r'-rL? 

4c7 'c U' (J jTI ) .LT.A“4^i‘ G0“ Tir”490 ’ 

4c= .“L“ ( ; 1 = 4^>, ( n *7 ls 
4 i ^ C . ^ T I Nu, 

-J -- I . 

r 

C dCA Ar.tA 'Xoj'USTNt NT F 0S“ H I'NG'TI P HAKE 



L, 


c 

c 


t '.fA'iUbi ^JL'JIN 


!fiA>': UT.crj .j 71 oLb' 

a:j-,T ,l4u TIP U* PJ-< hiNG ^OQT hake 

'll l~ [ A' I J J « I } . 4 T • A 4 A% I ltd Tu Ol2 
GC 73 jia 

517 Jj= 

IFCJJi o14*o14ij1o 

514 .4l;( 2 ,n = c 

If lni.HtI).^o.u) GJ TO t.16 
■rPIT!ia><tl55) 111 ,.4LWII I ,«LCI 2,1 ) 
bTO? 

■5lt CJVTl4T3t' 

530 CC4TIhU? 


IF IASnC.CT.O.C ) 
T?' 


Go TQ 700 


•KrTTre- 



InIIIALIZAT ICN'FUik WING ROOT KAXt 
K 2=0 

> B sO . 4 

YR=0;‘C. 

Re2=XTCC(1I 


PLHFH 

211 

PLHFH 

212 

PLHFH 

213 

PLHFH 

214 

PLHFH 

215 

PlHfh 

216' 

PLHFH 

217 

PLNFH 

218 

PLHFH 

219 

PLNFH 

220 

PLNFH 

221 

PLNFH 

222 

PLNFN • 

223 

PLNFH 

224 

PLNFH 

225 

PLNFH 

226 

PLNFH 

227 

PLNFH • 

. 228 

PLNFH 

229 

PLHFH 

230 


FLHFH ZTT 

PLHFM ■ 232 

PLMFH ' '2^3 

FLWFM 234 

PLHFH' ■ 'ZTiT 
PLWM 236 

"TTRFTi 2TT 

' PLHFH 238 

PLHFN 239 

PLNFH 240 

PLNFH 241 

PLHFH 242 

PLHFH 243 

PLHFH 244 

PLHFH 245 

PLHFH 246 

PLHFH 247 

PLHFH 248 


nrwTT 

2^ 

PLHFH 

250 

plnfm 

251 

PLHFH 

252 

PLHFH 

253 

PLHFH 

254 

" ■ ~ "FThTH" 

-25? 

PLHFH 

236 

PLHFH 

257 

PLHFH 

258 

PLHFH ■ 

259 

PLNFH 

260 

' - 

zsr 

PLHFM 

262 

PLHFH 

263 

PLHFH 

264 

PLHFH 

265 

PLNFH 

266 

PLNFH - 

ZVT 

PLNFH 

268 

PLHFH 

269 

PLHFH 

270 

- PLHFH 

- “271 

PLNFH 

272 

PLMFH ■ 

Z7T 

274 

PLHFH 

275 

PLNFH 

276 

■ ■ PCWH 

■ 277 

PLHFH 

276 

PLHFH “ 

279 

PLNFH 

280 



YE^,«Y^0C^1) PDiFM ?gr 

X2»XEi ' PLHfM 282 

„ . -pCUFff TSJ 

*SHG=-1.0 ■ , PLMFM. 284 


c 

• -• 

JLT»i ; :: 

CO TO 220 

PLNFH 

PLNFH 

T¥F 

286 

ERROR MESSAGE 

PLMFN 

287 


630 

K«450 

PLNFH 

288 



A«C 

PLNFH 

289 



"2-Z2 

PLNFH 

290 



WRiTE(Irf,15C) K,2,A 

PLNFH 

291 



STOP 

PLNFH 

292 

c 



PLNFH 

293 


700 

CONTINUE 

PLNFH 

294 



IFIJRITE.EO.OI GO TO 760 

PLHFH 

295 



HRXT£(1H,65J 

PLNFH 

296 



tF(NEN.EO.l) hXlTE ( IW,b6) 

PLHFH 

297 



IF(NEH.EQ.2} hRITEUN.67) 

PLHFH 

298 



DO 750 1=1, L ■ . ■ 

PLHFH , 

299 



JL^HLUl 

PLHFH 

300 



K»4L/5 

PLNFH 

301 



IFl JL.Lt.5#K ) 60 TO 710 

PLHFH 

302 



K«K + 1 

PLHFH 

303 


710 

HRIT£IIH,70) I 

PLNFH 

304 



WRITEIIH,72) FLC«2,I> 

PLNFH 

305 



IF (HLHm) 7 iS,725,720 . - 

PLHFH 

306 


715 

MRITE (1R,75) FLKU) 

PLHFH ■ 

307 



GO TO 725 

PLNFH 

308 


720 

NR1TE<1>I,80> riHU) 

PLNFH 

309 


725 

DC 740 JJ=1,K 

PLNFH 

310 


740 

HRlTt f (J,A k( J, i) , J=Jj, JL,K) 

PLNFH 

311 


750 

CCf>Tl6Ut 

PLNFH 

312 

c 



PLNFH s 

313 



JL=«LIL» 

PLNFH 

314 



itRITE (o t 1'HG ) 

PLNFH 

315 



RRITEIIf ,1*15) ( J.Xut 1J> , xf E ( J) , J=^l, JL ) 

PLNFH 

316 

c 



PLHFH 

317 



SUf.=2.^SUN*G2 

PLNFH 

318 




PLNFH 

319 


760 

StETORN 

PLHFH 

320 


65 

FCR*'ATtlHi,2GX,47HNON-lIhENSI0NAL flCX-AREA 0 ISTRI BUTl ON ) 

PLHFH 

321 


66 

FORMAT! IH+,57X,21H - ( PHYSICAL HIHGJ//I ' " PLNFH ‘ ' 

322 


67 

FCRMAH lH*,57X,tlrt - (TRANSFORMED NIMG)//) 

PLHFH 

323 


70 

fokmaT(//2x,I2,ivh-th spanhise column/) 

PLHFH 

324 


72 

FORMAT! 12X, i6EF IkST i_.c. BOX =-,I3»22H-TH BOX FROM NIHC ROOT) 

PLHFH 

325 


75 

F0RMATU2X, loFf IR5T HAXE BOX «,I3,22H-TH BOX FROM WING ROOT/) 

PLNFH 

326 


80 

FQRMAT!l2x,ifcFLAST MAKE BOX =,I3,22H-TH BOX FROM MING ROOT/) 

PLHFH 

327 


85 

fCRMAT15X,5( lX,13»lPcll.4n 

PLNFH 

328 


90 

FORMAT!//1Cx,25mH!N(, iREA CALCULATED FRQH/ 

PLNFH 

329 



15X,2tHLEADING AND TRAILING EDGES -.lPEll.4/ 

PLHFH 

330 


< 

15X,26HSJnMATlON OF AREA OF BOXES -,1PE11.4» 

PLNFH 

331 


96 

fORNAT(//10X,43HNEGATIVE BOX AREA EXCEEDS ALLOMABLE LIMIT (« 

PLH^H 

332 


1 

1PE11.4,1H)/10X.3HAR(,I3,1H»,I3,1H,,6HHEM» >«1PE11.4< 

PLHFH 

333 


2 

5X,21HCCMPUTATI0N CONTINUES//) 

PLNFH 

334 


140 

F0RHATtiHl,9X,54HLEA0Il!G AND TRAILING EDGE POINTS AT EACH CHORDKIS 

PLHFH 

335 


SE ROH/13X,lHJ»fcX,7HLEAQING,8Xi8HTRAlLING/l 

plMfh - 

336 


145 

FORMAT!10X,I3,iP2£15.5) 

PLNFH 

337 


150 

FORMAT! 1H0,4X,29HPLNFM — NEGATIVE VALUE NEAR SNv I4«2X*A2«1PE11.4) 

PLHFH 

336 


155 

FOR«ATUHO,4X,25HORDER OF MING TIP BOX IN . 1 2,43H-TN- SPANHISE C(H.U 

PLNFH 

339' 


IHH IS HOT PROPERLY OEFIHEO// 

PLNFH 

340 


2 

15X«8HHLH(HEH,,12«3H) ••13/15X,8HMLL<H£M»tIZ«3HI -.ISt 

PLHFH 

341 



END 

PLHFH 

342 



ORIGINAL EAGE 
OF POOR O.UALTTV 




SUBROUTINE PO72(H2fHO«l1O»CK«0*A) ~ . 

ran 

•Z 

c 


- 


3 

c 


Tht VELOCITY FIELD OF A UNIFORM DOUBLET DISTRIBUTION 

'4ctfT2 

A 

c 


OVEK A BOX IS COMPUTED AT ALL POINTS AT HHICH IT HILL BE 

R0T2 

9 

c 


NEEDED AND SICREO IN THE ARRAY A IN CONHOH 

P0T2 

6 

c 



R0T2 

7 

c 


M0«N0 CONTROL THE NUMBER OF VALUES COHPUTEO 

P0T2 

8 

c 



P0T2 

9 

c 


M2 IS THE RAHCE OF THE SECOND SUBSCRIPT IN THE ARRAY, 

P0T2 

10 

c 


DImEnSIONLO Afz,M2,N2). BUT TREATED HERE AS AN ARRAY 

P0T2 

11 

c 


with tbC subscripts 

P0T2 

12 

c 



P0T2 

13 



dimension A(2,1) 

P0T2 

14 



H-nO 

P0T2 

15 



S*NO 

P0T2 

16 



OK=CKAO 

P0T2 

17 



Dk2=«0K*»2 

P0I2 

18 



Ni=R-l . 

P0T2 

19 



DK8=3k2/8.0 

P0T2 

20 



DK-<i»2.0*DKo 


2l 



DA1; = Ok.2/12.0 

P0T2 

22 



CM=0.5 

P0T2 

23 



OH=Jk»0 .3 

P0T2 

24 



Oh=O.S*Dri 

P0T2 

25 



CC=2.C»OR 

P0T2 

26 



CEm=OD 

P0T2 

27 



C1=0.23*0K2 

P0T2 

28 



tJS=CK2/24.0 

P0T2 

29 



DO i 1=1, M 

P0T2 

30 



&1=0.0 

P0T2 

31 



= 2 *0/ Ort 

P0T2 

32 



82=ji/34-un 

P0T2 

33 



23=-0.3*85 

P0T2 

34 



u3=3H*o4 + £!> 

P"0T2 

35 




P0T2 

36 



DCA^Z.G^OA 

PCT2 

37 



CN=1.C 

P0T2 

38 



K = 1 

P0T2 

39 



C3=0.0 

P0T2 




C4=C.O 

P0T2 

41' 



C7-0.0 

P0T2 

42 



CS=6.0 

P0T2 

43 



00 2 3=i,N 

P0T2 

44 



Al-OM/CN 


45 



C1*CM^ COSIAl) 

P0T2 

46 



C2=-CM* SIN(Al) 

P0T2 

”47 



C5=CM*CIN«Ai,C6) 

P0T2 

48 



C6*-CM*Co 

P0T2 

49 



C9=*Cl-C3 

P0T2 




C10»C2-C4 





C11=C5-C7 

P0T2 

52 



C12=C6-Ca 

P0T2 

53 



Atl.K|»B3*C9-e4*CL0-B5*C3-Bl*Cll-B2*C12 

P0T2 

54 



A«2,KJ=E-S4C9,e3*ClO-B54C4+B2«Cll-ei*CX2 

P0T2 

55 


23 

C3«CX 

P0T2 

56 



CA-C2 


tr 



C7=C5 

P0T2 

58 



C8>C6 

P0T2 

59 

« 


81-Sl-Dl -- 

■P0T2 

60 



B3-33-03 

P0T2 

61 

r 


'■ - 

P0T2 

62 



0A>0A«DD4 

P0T2 

63 



CS^tf*?75 

P0T2 

6*4 


2 


P0T2 

65 



cr*-tM*i.o 

PDTr 

■' ■ 66 



OM^OM+DOM 

P0T2 

67 


3 

OfiM’ffDHTDD 

P0T2 

68 



no ■; Il-1»2 

P0T2 

69 



Al-l 

P0T2 

70 


hy 





00 5 J«1.N 

P0T2 

.71 


“Off" 4 T= 1 ,751 ' r - - 

FOT? 

■ 72 


K'Kl+t'-l 

P0T2 

73 

4 

A< lL»K)>AaL.K)-An'L»K-n 

P0T2 

74 


AUL«Ki)=2.0*A<lL.Kl) 

P0T2 

75 

5 

Kl*ki*M2 

P0T2 

76 


cn»o.o 

P0T2 

77 

“ “ 

“OH '.OVD 

~P0T2T 

■ 7ET 


00h=DK 

PDT2 

79 


DC 12 I =1,11 

P0T2 

80 


C7*0.0 

P0T2 

81 


C8«0.0 

P0T2 

82 


C9»0,0 

P0T2 

83 


CI0=0.0 

“POTZ 

84 


P1«0.0 

P0T2 

85 


PZ«=0.0 

F0T2 

■ 86 


CN-1.0 

P0T2 

87 ■ 


ff6=0V5*D!CK ' _. .. - -. -- — 

Ptrrz 

80 


K«I 

P0T2 

89 


UU 10 

PU 12 

—90' 


Al-CH/CN 

P0T2 

91 


A2=DN7CN ~ •' ■ ■ “ 

nrrz- 

92 


IF tAl-0.21 7,7,8 

P0T2 

93 

7 

0I»Z.D-Ar?*I73.O' ” 

~ T>irT2“ 

94 


B2— 0K/«6.0*CN) 

P0T2 

95 


GU TU“9 ; 

PDT2 

9b 

a 

B3« SlHtAl)/Al 

P0T2 

97 


B1»2,0»B3 ^ 

ponr 

9S~ 


e2“(83- COSJA1H/A2-OH/CN4B3 

P0T2 

99 

9 





84« SIN(A2)/CN 

P0T2 

101 


C3-»i«B3*B2404 

F0T2 

1U2 


C4-82*B3-B1*B4 

PQT2 

103 


■^■?-OH*CN 


■HHIIIIMtLV 

1 

C1=B5*C4-2.0*C3 


105 


C2»-2.15*C4-B&*C3 ; 

P0T2 

to* 

1 

C5-C1-C7 

P0T2 

m 

j 


P0T2~ 

' aM 


P3-P2-864CN 

WT2 

1%9 

S 

P4»f 3*2.040K124CCH-1.0I 

Ifl'Jt] tM 

— =TP|Sr 


A C 1 ,K > » A < 1 ,KI 14C6H*3*C3~P44C9 

►012. 



A(2,Ki«A<2,K}«Cb«Pl«C5*P3»C4-P4«C10 

— P0T2T 



P1-P1-H)H 

**0?2- 





9S^S9JLm 


CN=»CN*2.0 




C7-C1 


I 



C«*C2 

C9»C3 

CJ0»C4 

.Bfc»B6*8Ki2 

TolRfHtz ■ 

Sttsim^PH 
12 ODn^DDHtDD 








->ji 


V;-' " 

' '»\r,*. 


03*C</ ( 2 .«4 3.' 141 $92451 
M-H2HI 




K«1 

Ai>0.0 


-■If": 


00 14 1>L*N 
Cl-03* S1WI41I 
C2-- 034 COSIAl) 
DO 13 1-1»W 


• , . ,; ' 


DFE -A<l»K)4ei»AC2<X1*C2 
A1 2»KI»Al2.m»Cl-All»XI* Cg 
AU»K»«DFE ^ 

13 K»X-»i 


K>K>M1 
14 A1>A1-»0H 


RETURN 

END 


^8 







n, n r> o rt 





SUSTtaL-riNt CREDt Sf O/.SFDT » SlDh, KSFD, SFMX»SI-l'V,SF«H,KSFn,DA»T 

DREO 

2 

I 

. ■ ,XJ<.YY>I>nL*HLC*»l»hu.NE,NEti,NH,KSFS,.lP(UNT» 

OREO 

3 


r>l>i£'<SION SFOy(«B,U,SFDY(NB»nJ),SFOriCN«,ML<l.KSFO<l} 

DREO : 

•• 4 • 

I 

,5t-.4XaJ»5f.iY<l),SFt<h(l).UA(l) , T«NM,1 ) »rtLC(2tU 

OREO 

5 


,/.l 1 1 J.YYl 1) 

OKED ' - 

6 



OREO 

7 


£RKR-l.t-C£, 

DREO 

8 


ieiffc*IP«iNT/lC 

OREO 

9 


jMic=I?RI«T/IGO 

DRED 

10 


IFIMtH.LO.il IRlTt » IPR IN T-10» I R I Tl 

DREO 

11 


4SFC(MM0A<Sjs ) 

DREO 

12 



OREO 

13 


IF INP) 73C»17C»100 

OREO 

14 

■ 


OREO 

15 


A SPLINt-SURFACt Fjk the DtFLtCTION IS FITTED TO VALUES 

OREO 

16 


OF deflection at OIVt«* POINTS. 

OREO 

17 



DRED 

18 

IOC 

IF <NS-NP) 73C,120,123 

DRED 

19 

120 

CGF'TIM.E 

OREO 

20 

~ ~ 

iP» iOO 

' OREO 

21 


OC 14C IP=i,NF 

DREO 

22 


SFCXC IP ,n)*ilA<KP + i)/DA(^4) 

OREO 

23 


SFCYl tP»«) = 0A(^P♦2)/JA(2'tJ 

OREO 

24 


SFOrtI IP»«>=uA (KP*31 

OREO 

25 

140 

<P»KP«4 

DREO 

26 


SPLlNE-SUKFACt FIT DATA 

DRFB 

TT~ 



OREO 

28 

145 

CCnTINUJ ' 

ORED 

29 


IFIIRITE) 150,lo0,i50 

DREO 

30 ■ 

150 

HR I ft 1 IH,30T~F 

OREO ■ ■ ■ 

■ '31' 


KFJTEUMtiO) 

DREO 

32 ■ 



OREO 

33 ■ 

loC 

CChTINOE 

OREO 

34 


CALL SuyFlTFn ,KS?D'(rt).r,STOx(l.fi) ^TFCYTl ,M),SFDff(I,ffI.lRITEJ 

■ OREO 

35 


I.C TO 200 

ORED 

36 


PRESENTLY FOR PHTh AND PLUNGE OT'FORH Z-AO*Al»X 

“ DRED 

■ 37 

170 

SFiHt l,m«DAJ52) 

ORED 

36 


■ S‘F L rTr2TS' J -DA I 5 3 WAT 2A » 

DREO 

■’39 


SFCN(3«H)=0.0 

DRED 

40 


GO T3 14? 

■ DRED 

41 


200 ICh£CX=C 


ORED 

42 

RETURN 
730 IPR-98 


"DRED 

OREO 

44 

■ GG TO 750 


DRED 


750 WRITE 


OREO 

46 

C (Tw.aOllPR 


DRED 

ajT- 

STCP 


OREO 

48 

io FCRIu7llHO»10X»t3HCO«PUTEO DEFLECTION - 

A04Al*X*A24y*“ SOM OF 

■RTTT*nDRED 

vr 

HRa»**2)*ULCC<Rm*421l> 


OREO 

50 

’2o fDR*iAmWfl;iOX,l6H&ft£D — 6A& 0ATA,15» 


DRED 


30 FORNATIlHOi8X«26HPriVSICAL PLANE — MODE 

HO. *13) 

DREO 

SR 

END 


ORED 

53T 


43 



$U8R0UTIM£ HV||,4^«Vf,X33^yiKmttHt^SFOX«SFI)¥*S»m«K$FD«S«0»at 

CALCUCATFS OOXMlfASFI VCtO&ftV OFISimSVTlQM tR£«& AHO IMAClIUCyt 
01 HFMSIOK XX t AikXV I itVKLl II 
* .SFeX4HaslF»$F{>VlfiB«M>l9SF^fXH*tt0i*KSF0tl» 

“TSTK^^imSTTir : — 

lFtlt£8.£<l.ll 





4^ 


IT 




— g- 

IFtKSF0iKl).E0.O.AttO.l(EH.aa«Zl JLF-KSFS 

HVAL 

10 


~9f^ 

11 

00 5 

WVXL 

12 

00 5 i-l,JL ^ ■ 


VS~ 

00 5 K»l*2 

»VAL 

14 



Hr 

00 80 I.1,L 

MVU. 

16 

3C=HtTn ~ ~ 

Bwo:: 

nr 

IMiLl 40tOO«40 

HVAL 

18 

40 CONTINUE 

BVAi: 

pr 

DO 70 J«1,JL 

HVAL 

20 

It* 1 XX f i ) *0 K « IXlfcl ji 'tiD W (>U tU 7Q 

w?*c — 

ZT" 

CALL SURF2(XXIi>*YYti>tl*lfl«VALU«VALU0,0HySF0X(l,Nlt 

HVAL 

22 

» nrFBY(l*Nl>tSFt>Hll«ai),JLP«2l 

“ HVSt 

23 

S(1«J.I1«VALUC 

HVAL 

24 

SC2«JtI »>CK4 vALU 

BV*t 

25 

70 CONTINUE 

HVAL 

26 


HVAL 

ZT" 

HVAL 

28 

HVAL 

— 25 

HVAL 

30 

:QTE;' THAST'AHtCETr'WAl. ■ 

- . 

3/> HVAL 

32 

HVAL 

■' “33* 

HVAL 

34 

HVAL 

35 

HVAL 

36 

HVAL 

37 

WVAL 

38 

HVAL 

35 

HVAL 

40 

HVAl 

41 

HVAL 

42 

HVAL 

43 

HVAL 

44 

HVAL 

45 

HVAL 

46 

HVAL 

47 

KVAL 

■ 48 

HVAL 

49 

HVAL 

50 

HVAL 

. 51 

HVAL 

52- 

HVAL 

53 

HVAL 

■ 54 

■HVAL 

55 

WVAL 

56 

1) .S<2,JltI),S3.S4 HVAL 

57 

HVAL 

56 

HVAL 

55 

HVAL 

60 

HVAL 

61 

HVAL 

62 

HVAL 

63 

HVAL 

64 

HVAL 

65 




12 


IFIIKITEI 100«200»100 

"tnOTE (IHVlff} 

MRITEI1M*12) Ml 


120 

'ITO 


F0RHAT(lri*,56X,28H PHYSICAL PLANE HQOE N0,,J 

DOTT'O I-l.L' 

JL="L {1 ) 

■ “ir“m.TTr(5TnrruTi — - - 

110 HRITi- n«*,20) I 

JLP=JL/2 

IFC JL-2SJLP.NE.01 JLP=4LP+l 

CO l&O J«1 ,jlF — 

Sl^S3RTiSJl,J,I)SSll,J,lJ+S£2tJ f nss(2« Jtin 
IFtSl'.CT'.G.Or'CO To 12 ‘J 
S2»C.O 
GO TO 13'T 

S2=57.29S73SATjn2{S(2,j,i ) ,S U , j,I J) 

ji=j»Jlp" 

I F ( .LL . Jl i CC TO l30 
Jl«0 ■ 

S3=0.0 _ 

S4=C'.0 
GO T3 ItjO 

■nriiTTKUE 

S3-SORTiSJ_W J 1, I J_*_S li..Jl.Il+S(2,Jl,l)*S(2,Jl,I>) 
1 F i S 3 . C T . C . 0 ) Go T u 
S9=0.C 
GC‘ to' 16G 

S^ = S7.29S7 b»ATAN2< j ( 2tJitl),s'(l,Jl,l)) • 

•<R1 T£« I(i, 23 )J«S(x,Ji: ),S'I 2,J,I ) ,SltS2tJl,S(ltJl, 
CJDN T I NUE 

THESE ARE THE LPMASriES 


T5^- 


IS5 

‘160 

170 


200 CONTINUE 


RETURN 


20 


25 


FORHATCIH »5X,2{2Xtl3, IP9E13.5J) 
END 


50 



oo’oH'o'J J noonn 


ORIGINAL Bags jg* 
OF POOR aUiBITV 


SLBkJl,Tl;St .-<C)iF(XX,Yr ,XTD 0 fY 10 G,XTE,XLE «AyAR«0AN*TtEOG«S9CK*0*YE 

r ■ ■ ■ ■ ~ i‘i ('TInt »MLT,rtrc » am* : ^ 

2 ,JM 4 x,t.nL,MLH,hl,NtM»NST»lW»Me,MC 9 HE»M 0 i' , . _ ' 

ICMtOX^O - P,L7thTl4L AS CQMPUTcfi'lh 80XPC ' ' • ■' 

lChtCn = l -W!Ih LlACIni, LDCE CO ftREC TlOM ■ 

IChcCK=<’ - f-^issukt S'JtFf IcIEMT COMP.uTeO from corrected POTEHflAL' 
ICHtCK-3 - PKlSSURL CJEFEIUEMT COMPUTED FROM TRANSFORMED VALUE . 

For physical hing - ~ 

D1 MENS ID'S XX ( 1) ,YY( I) , X T Oo < 1 ) 9 Y lOCU > i X TE ( 1 ) « XLE ( 1 ) i A (2 iMCtHB I 
1 . AL ,^-OG t 1 1 ,S{2»nB,MBI ,T IME-.M-Eiil ' . ' 

■2 «nUl)>r!L»<m«MLT(29i) «rlLC(2«HBi«AMA(HB>Mei ' 

IPk 1N=IPRINT 
KRITc = IPKlri/iCCGo 
IPRIN = IPR1N-1 CCO()»KkTT£ ■■ 

JRITl =1PKIN/1C0 

lPRIN=lPRlN-lCC*Jkl Ec ■ - • 

IRITc=IPRlN 

ifinlw.ll.i) Cl to I 1 O ' ■ ■ 

<«lTE=KRnE /IC 

■ jm-m j^rTTCVor — 

IRITE = I!UTi./lC 
GC T(, 150 

lAO KRlT>:=i<k.nt'(KRIT6/lO)*10 

jRife=jRncG( jkrTc/u»rpr6' 

IRITE»IkITt--( IRlTE/lu)PiO , - ' - 

■'l30 CCNirNUc ' “ ~ 

NS«0S»0 

DHiC,5*J ■ ' ' ~ r ■ 

CAlL~ e 0 X pa { X X, V Y ',TT0 G ,YT0G»XT£»A,AR» D an VT r57CK7tr;'AMS'" 

_ 1 1 k»i-,rtmhCiM L< MLC«MLH»Hl9NEV<iHST9 JRITEI 

eOXPO CCMPUIhS THE POTENTIAL VALUES IH EACH BOX. 

"THE^Y ARE~ STOjivtC iH "Tr*!; Ai^RAY' S* 


BOXP 
"BOXP 
BOXP 
'BOXP ■■ 
BOXP 
BOXP 
BOXP 
BOXP 
BOXP 
BOXP 
BOXP 

“BTfif? 

BOXP 

eoxp 

BOXP 

B'OXP 


BOXP 

20 

BDXP 

21 

BOXP 

22 

■ BOXP 

" 23 

BOXP 

24 

BOXP 

“25 

BOXP 

26 

■ BOxP"' ■ 

27 

BOXP 

28 


ICHEtK=0" ■ ... 

MEH=MLM 

JflXK^'JFrO" 

IFllRITt.EO.o J GO TO 270 

HTUTl (ITit 30 ') ■ ■ ' ' ‘ “ 

200 DO 205 I=1»L 

DO 205 J=lt‘jMAk ~7 

DO 205 X=1.2 


I iJllI 

215 MR1TEUH,32) Ml 

30 FORNATUHl. rcXjfalHPOfENHAL'CATCULAtED' (REAL* IMAGIP 
4 PHASE ANGLE)) 

“ 32 ~F0RMAT ( IH* f IX ,28ri---rri‘YrrCAL' PL ANT— -MODE NQ.* 13/ ) ' 
C PRINT-OUT 


218 DO 250 1 = 1, L 
JL=MHI) 

. ' “ IF <3L) '220,25C,220 
220 HRITE(IH,20) I 

_ IF(JL-2»JLP.NE.0) JLP=JLP*1 

00 250 J=1,JLP 

S1»S0 RT(T ( J , l, l) » TtJ,l, l)*T(.) , ) » T('j, 1,2) ) . 

if (Sl.NE.O.C ) GO Tl 225 ' ■ . 

'S2=0.0"' 

GC TG 226 

225 CONTINUE 

$2 = 5 7.29578 A AT An 2( TlJ,I,2l,T(j;i,D) 

226 Co'nTi'nuE 

Jlr-tL.-lLl’ ; _ 

If <Ji .LE. JL) cc'"lt. 

J1=0 

53=0.0 

,1 1 =e.c 


■ BDXP 

BOXP 

■-fflWP 
BOXP 
'■ BDCT 
BOXP 
BDXP ~ 
BOXP 

■ BOXP 
BOXP 

•“5DXP — 
BOXP 
'BOXF 
BOXP 
•~BffXP — 


BOXP 

45 


Jjy- 

BOXP 

46 

; "BOXP~ 

~\7 

BOXP 

4B 

“■ BOXP 

49 

BOXP 

50 


5T 

BOXP 

52 

BOXP 

53 

BOXP 

54 

BOXP 

55 

BOXP 

56 

BOXP 

. 

■ BOXP 


8QXP" . 

'59 

BOXP 

60 

BOXP 

6l‘ 

BOXP 

62 

BOXP 

63 

BOXP 

65 

BOXP 

65 

EOXP 

66 

BOXP 

67 

BOXP 

68 

'BOXP 

69 

BOXP 

70 


51 





TIJl , I , i )-C.O 

BOXP 

■71 



CO TO 240 . . 

BDXP 

■ .72- 


230 

tONTllUE • 

BOXP 

73 




BOXP 

74 



CC rj 234 

BOXP ■ 

75 



S^=G.O • ■ . 

BOXP 

76 



GO TO 236 

BOXP 

77 


“F34 

CCNTIhUt 

■ BOXP 

78 



5<. = t V - 1 T ( Jl, I ,2 ) » I( Jl , 1 ,11 1 

BOXP 

79 




BOXP 

80 


240 

4RUM1",2-. ) ^,nj,I,ll,T(J,I,2),51,S2 

BOXP 

81 




BOXP 

82 


250 

COM I.4Ur 

BOXP 

83 

c 



BOXP 

84 


27C 

CCMINUE 

BOXP 

85 



If ( lCt-£CK.Gt. 11 Go TU 400 

BOXP 

86 

c 


LEA0I.5C iUGfc CLR%coIION 

BOXP 

■ 87 



OQ 3t> 0 

BOXP 

88 



IF INLT < 1, J 1 .EC.O) Ui 285 

BOXP 

89 



Il=rLTli,J) 

BOXP 

90 



I2=nLT(2, Jl 

BOXP 

91 



S3“0.0 

BOXP' 

92 



00 280 1=11,12 

BOXP 

93 



K=1 

BOXP 

94 



IF (A»14,I1.G£.1.01 Gd TO 290 

BOXP 

95 


280 

CONTINUE 

BOXP 

96 


285 

S3-^C 

BOXP 

97 


290 

IFlK.ea.U S3 = C 

BOXP 

98 



S3=XXlK)-XLt( J)-0ri-»S3 

BOXP 

.99 



OG 370 1=1, L 

BOXP 

lOD 



Sl-lXX (I l-XLE ( J1 1/53 

BOXP 

101 


■» 

IflSl) 30C,34C,3i0 

BOXP 

IBT' 


300 

S1=0.0 

BOXP 

103 



GC TO 340 

BOXP 

104 


310 

IFITE) 320,32C,330 

BOXP 

105 


320 

Sl-S1*(XX{I)*XLE<J» )/(S3*2.4XLEU)> 

BOXP 

106 


330 

IF! Sl.GT.l.C) 51=1-0 

BOXP 

107 



Sl-EDCI J»*SCRT1S1) 

&0*P 

ni5~ 


340 

DO 350 K«l,2 

BOXP 

109 


350 

SlK.J.t 1 = S14S(K, J,1 1 

BOXP 

1X0 


370 

CONTINUE 

BOXP 

111 


380 

CONTINUE 

BOXP 

xii 



ICHECK=-1 

BOXP 

1X3 



IFURITEI 390,400,390 

Bbx4» 

iPT 

- 

390 

WRITE cm, 36) 

BOXP 

U9 


36 

FQRNATC1H1,10X,61HPOTENTIAL corrected (REAL, IHACINART, ABSCMLUTE, 

BOXP 

• XU) 


S PHASE ANCLE )l 

B02F“^ 


c 


TRANSFER TO THE PRiNT-OUT OF CORRECTED VELOCITY POTEHTIAL 

BOS> 

Xt« 



CO TO 200 


~~~~SW 


400 

CONTINUE 

BOXP 

■ -12» 



IFIKRITE.EO.OI CO TO 700 

80XP 

IBS' 



IF(lCHeCK.&£.2) 60 TO 520 

BOXP 

122- 

c 


CALCULATE AND PRINT PRESSURE COEFFICIENT 

BOXP" 

■' • 



ICHECK=2 

BOXP 




DO 410 I-1,L 

Box# 

i25 



JL<=NL(I) ' ' ~ 

BOXP- 

126^ 



00 410 J-l,JL‘ 

80XK 

IzS^ 



T(4,l ,3)=S(1,J,1) ■ - ' 

BOXP 

120 


410 

T(J.I,4)-S(2,J,I) 

BOX# 



420 

JNAK= JHAX 

BOXP ' 

130 



00 430 I»1,L • ' 

■^OXP" 

> in 



jl<=hl(I I 

BOXP 

132 



00 430 J>1,JL 

BOXP 

or* 



DO 430 K>1,2 

BOXP 

134 


430 

T< J,2 ,K)=0.0 ~ = — ' — 

BOXP 

T35~ 



00 500 J«1,JHAK 

BOXP 

130 



11>NLT(1,J1 

■B0XP~ 

137 



I2-»-Tt2»J)~NLTCl,J»+l 

BOXP 

138 



imn 

"BOXP 

I57~ 

' 


1K«12*1 

BOXP 

140 


52 



ORIGINAL BAG'S m 
OP POOR auAiiiTy 


. c 

TTl, 6,Sj»xLEUJ ^ ^ ^ = — 

TM. 8*51»XLtU» 

“BOTP 

BOXP 

— KX" 
142 

m.i0T3i»-CiT0 ■ ' ~ 

Til,12,5)-0.0. , . 

BOXP ■- 

143 
' 144 

If- ' 1 ABS IXLt f Ji-XXt 1 i ) ) .GT« !• E— 05) GO TO AAb 
- ADJUSTMENT ■- LEADING EDGE ANO FIRST BOX COINCIDES ' 

-BIIXP 

BOXP 

*145“ 

146 


^ 


147 


IK- IK-1 

BOXP 

148 


A«'CBRTTNOE ” “ • — 

BOXP 

"T4^" 


DO 450 l-2« IK 

BOXP 

150 


— m;“6 .sr^“rm - — ; - . - - 

BT5XP " 

■ 151 


TII, 8,5I«XXU1) 

BOXP 

152 


T ll Jvll *3 J 

■^BXP” 

153 


- T(IVi2»51-TI J,II ,41 

BOXP 

154 


— ni-TTirsTiofxin) ■ 

BOXP 

155 


450 n-ii-1 

BOXP 

156 


CALL SF'LNiU 2',T't l'» 1 ,'5Tl' Vt'l 1 » 2 » 5) »T ( 1 , ' 4, 5 1 

BOXP 

157 


% »IK,T(1. 6<5)yT(ltI0t5)«TU,14*5')«NSM0S) 

BOXP 

158 


OOTTFCnTT r2TTTlT~T75TrnT» ^.'5)»TTTT 5,TT ' , ’ 

BOXP 

■IBO 


» ,IK,TU, 8«5) «T(1,12«'5) iT<ltl4f 5)«NSMaS) 

BOXP 

16 G 


u«rr ' ■■ 

BOXP 

Ifal 


OC 470 1=1,12 

BOXP 

162 


T { J ; 1 1 , 1 > =f r I V4' , 5"l -CK*T (1,3,5) 

BOXP 

163 


T(J,II,2)=T(I ,5,5(+Ci<«HI,2,5) 

BOXP 

16 4 


4 7o~n i rm - . . 

BOXP 

165 


500 CCNTIHUt 

BOXP 

166 


HRltc( IH,40)‘ 

BOXP 

167 


40 FORMATI IHi, 10 k,61HPkESSURE COEFFICIENT !i.c*L, IHACINARy, ABSOLUTE, 

BOXP 

168 


1 PHASE ANGLE)) 

BOXP 

169 

c 

transfer TO PRINT-OUT SECTION 

BOXP 

170 


GO TO 215 

BOXP 

171 


520 CONTINUE 

BOXP 

172 


700 CONTINUE 

BOXP 

‘ 173 


RETURN 

BOXP 

174 


20 F0RHAT(1H0,5X,I2,19H-TH SPANHISE COLUMN) 

BOXP 

175 

1 

25 FORHATUH ,5)(,2(2X,I3,1P4£13.5)) 

BOXP 

176 


END 

BOXP 

177 


53 



r> o' o n o oo 


SUSKOUTIHE BO>i>0(XX«rr«XTDG« YTOCtXTEiAf ARtDAN«TtSfCK«OvAHA 

V rnr,'l.,HB«HC.HLtHLC*HI.U,Hi,HCU,HST7Jftim 

SULUllflH OF SIMULTANEOUS EQUATIONS FOR THE POTENTIAL 


SOXPO 

■BOTFir 

BOXPO 


2 

TT 



OlnEMSlOM XX( U«Yr( U,XTOC(l)TTTDG(I)«XfE(l) tA(2«NC*ll 

BOXPO 

5 


1 . ,ARtMB,H6l*r(2«HBfl>*S (2»HBf l>«ANAfMB«HB) 

BOXPO 

6 


2 ',NL<l)>HLCi2f'li«HLHIl) 

BOXPO 

7 


Ch«0.5*0 

BOXPO 

a 


bC«2.O*0 

BOXPO 

9 


Ll^^-2*nLlLJ 

BOXPO 

10 


IF I JRITE.EO.O) GO TO 25 

— B'oxptr - 

n 


lF{NEW.£Q.2.flt,«l.GT,l» CO TO 25 

BOXPO 

12 


print INFLUENCE CUEFFICIENT 

BOXPO 

13 


MRITE(IM>100) 

BOXPO 

14 


jl.2*Hl(l) 

Btrevo 



K-JL/4 

BOXPO 

16 


If ( JL-4<fK.NE.C) ii-K*l 

^TJXPD' " 

T7 


OC 20 I«1»L 

BOXPO 

IB 


1 1=1-1 

'Bflxnr • 

■ 19 


MRlTfcdH.llCI 11 

BOXPO 

20 


"DTr ,« ■ — 

BUXPO 

— 71 

15 

tfRITE IIH,12C) (( J.A(l,J,l)tA(2« J*in« J>J1«JL«KI 

BOXPO 

22 

20 

CCKTINUE 

BOXPO 

23 

25 

CChllNUE 

BOXPO 

24 


il=NST 

BOXPO 

25 


Artl=1.0 

BOXPO 

26 


CO 26 1 = 1711 ----- — 

■ “BOXPO 

27 

2fc 

ANtsANINKXTOGd 1 >AMI) 

BOXPO 

28 



BOXPO 

29 


1F( JRITE.EO.O) CO TO 2d 

BOXPO 

30 


hKIT£UH,l68j ■ ' ■ 

¥OXPO 

31 

28 

CChTiNUE 

BOXPO 

32. 


1 1^0 

~ ■ BOXPO 

‘ 33 


NFLKS = DAN 

BOXPO 

34 


QC •JO I = 1«L 

&OXPO 

35 


X=AXt 1 ) 

BOXPO 

3b 


ACJUST UPSTREAM INFLUENCE 

BOXPO 

37 


XO * 1 

BOXPO 

3B 


LH^ ■ bQ« U 1 bU [U 30 

BOXPO 

“39 


KO = «AXO<l»I-MFLNS*il 

BOXPO 

40 

30 

CONTINUE ■ ' ' ■ ■ 

■ BOXPO 

41 


JL=HL1I) 

BOXPO 

42 


IF IJL.EO.OI CC TO'OO ■ 

BOXPO 

43 


DEFINE wing and WAKE BOXES 

BOXPO 

44 


CF HA'KE BOXES IH 

mnr 'boxfo" 

45 


JE=N0. OF MING BOXES 

BOXPO 

46 


■jS-ORBer 0"F“FT«^T wITiC‘ Box ■ BtiXPO 

47 


JH-OROER OF LAST WING 

BOX BOXPO 

48 


jSW-OROER Of~FlTST' WAXE 

BOX’ ¥OXPO 

“ 49 


JHH«OftOER OF LAST WAKE 

BOX BOXPO 

50 


JN = 0 

BOXPO 

“51“ 


IF(X.LE.ANI) go to 34 

BOXPO 

52 


J«1 • 

“BTOXPTT ' 

53 

32 

IFtJ.GT.JLl GC TO 34 

BOXPO 

54 


IF (X.LE.IXTE ( J)-»DHn GO TO 33 

BOXPO 

55' 


Jh= JW4l 

BOXPO 

56 

33 

j.j«l 

5UWD 

57 - 


GO TO 32 

BOXPO 

58 

34 

Jc = JL-JN ■ ' 

BOXPO 

59 


IF OILnlll) 3t*35«35 

BOXPO 

60 . 

35 

JS *JW*1 

BOXPO 

61 


JSN-1 

BOXPO 

62 


GC TO 37 

BOXPO 

63 

3e 

JS =1 

BOXPO 

64 


JSn=JE *1 

BOXPO 

65 

37 

JN =JS +J(T-1 

BOXPO 

66 


JN».= JSW + Ji*-! 

BOXPO 

67 


IF Ili.L'J.O) CG TO 50 

BOXPO 

68 


SLST-^ACTliiN OF CoNTRldUT IONS OF PRECEDING ROWS TO UPHASh 

BOXPO 

69 


CG 47 J=JSf JN 

BOXPO 

70 
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DO 45 K=^0^ 11 

BOXPO 

71 




6QXP0 

72 



Kl^l-^l-H. 

BOXPO 

73 



IF f.L.EO.O) lD to 45 . 

BOXPO 

74 



DC -40 N*1,K' 

BOXPO. 

75 



IF INFH.EC.4) go to 3b 

BOXPO 

76 



N1«N-»J 

BOXPO 

77 



N2-I46S IN-JJ* 1 

BbXPO 

78 



A1 = AI 1 ,N1,|( I )-»*« 1»N2,K11 

BOXPO 

79 



A2*A(2,NltM )-*A(^«N2,Hl) 

BOXPO . 

80 



CO TO 39 

BOXPO 

81 


38 

COKTINUI 

BOXPO 

82 



TFHP=(AHA( J,I )«A«Alh,K) J/OD 

BOXPO 

83 



yHl=tYY(J>*YY«F) UTtNP 

BOXPO 

84 



VM2=(YY« J>-YY(M )*Tt«P 

BOXPO 

85 



YH2=A6SIY«2) 

BOXPO 

86 



N1=YH1 

BOXPO 

87 



Fh2«YK2 

BOXPO 

88 



YHl = YtU-FLCAT IM ) 

BOXPO 

89 



YM2-Yf12-fl-aAT (K2) 

BOXPO' 

90 



Kl=rNl + l 

BOXPO 

91 



N2»N2*1 

BOXPO 

92 



IF (Nl.CT.LIfl2 I GO TO 138 

BOXPO 

93 



4RN10-A(1,N1,K1) 

BOXPO 

94 



AIN10<A(2fNivKl) 

BOXPO 

95 



IF (N1*1.GT.L1H2) go TO 238 

BOXPO 

96 



ARNll-Atl,M-»I,Kl) 

BOXPO 

97 



AINll-=A(2,Nl4l,Kl> 

BOXPO 

98 



GO TO 338 

BOXPO 

99 


138 

ARN10-A(1,LIMZ,K1) 

BOXPO 

100 



AIN10=Af 2yLlH2,Kll 

BOXPO 

101 


238 

ARN11»A(1,lJM2,K1) 

BOXPO 

1 O 2 



AI«11=:A(2,LIM2,K1J 

BOXPO 

103 


3 3'8 

CONTINUE 

BOXPO 

104 



APNZO-AdtHZ.Kl) 

BOXPO 

105 



AIN20«A(2tN2iK1) 

BOXPO 

106“ 



ARN21=A(1«N241,K1) 

BOXPO 

107 



”AlN2i«AI2»H24l,Kl) 

BOXPO 

nre“ 



Al=ARH10»tARNll-ARN101*YMl*ARN20*tARN21-ARH20J*YH2 

BOXPO 

109 



A2>AlN104(AlNll-AIN10)*YHl-»AIM204{AIN21-AtN20l*yH2 

BOXPO 

110 



Al«A14AKA«Ni,K) 

BOXPO 

111 



A2>A24AHA(N,K) 

BOXPO 

112 


39 

CONTINUE 

BOXPO 

113 



HT-1.0 

BOXPO 

— rrr 



IF IN.GE.NLCl 2«K).AHD.NLCf 2tK).NE.O) HT>AR(NfKI 

BOXPO 

115 



SU«J*rt=SUi J*n-(Al*Sf l,N,KI-A2«S(2tN*Kl)4HT 

BOXPO 

116 


AT 

SI2<J«I )>S(2« J«I l-(A2«Sf l«N*i04AT*SI2«M*KI)«WT 

Bireps 

117 


45 

CONTINUE 

BOXPO 

u« 


47 

continue 

BOXPO 

110 

c 


SETTING UP KATRIX FOR SIHULTAMEOUS EQUATIONS 

BOl^ 

120 


50 

bO 55 

BOXPO 

KIT 



DO 52 K>1,JL 

BOXPO 

122 



“Ni-j+k 

BoxpO ' 

123“ 



N2-IAB5I 

BOXPO 

126 



HT-1.0 

BOXTO 

123“ 



IF (K.GE.HLC<'2tn.AHO.>1LCI2tI).NE.O) HT-AR(K«It - 

BOXPO 

126 

c 


\ • 

~BfiXPO — 

127“ 



IF( J.EQ.Kl MT-1.0 • . . - - 

BOXPO 

128 

■C 


- 

Boxpo 

129 



T(l>J«K)-(A(l,hl>I}»A(I,N2«ll>4HT 

BOXPO 

130 



-TI2»J»Kl-IAI2,HX,l)*AI2»N2»l>l*liT 

BOXPO 

131“ 


52 

CONTI NUE 

BOXPO 

132 

c 


SUBTRACTION OF tOHtRIBUT 1(W 

BOXPO 

TST 

c 


FROM HAKE BOXES — SI NINO- 

BOXPO 

134 

~c 


~ SIHiHGI-TIHAKEIAPHTIHAKEI 

Bmtpn — 

ISF 



IFIJH.EO.O GC TO GO 

BOXPO 

136 



DO 5b j=jS7jh 

““BOXPO — 

■ ITT 



DO 55 N=JSH,JNH 

BOXPO 

136 



~Sll*J>I)-S(l,JfIJ-Tll>J«N)4S(ltNyI|4TI2«J«N14SI2tNtI) ‘ 

“nsoxpo — 

T5V 


55 

S(2«J«I >-SI2»J,Il-Tll,J,NJ*SI2,N,I)-TI2»J»Nl4SIliN«I) 

BOXPO 

140 


55 




•>6 COHUNUE 

bOXpo 

— nrr 


IF im.Hm.LT.0) GO TO 60 ■ 

BOXPO 

142 


■ - ft6-P051Tn3N ILCTEHTSTOrT BUXWT ' 

- .x«- 


OQ 59 

BOXPO 

lAA 

' 

-NTT-irf jh ' — 

•BDXPB" 

«5 


00 58 

80XP0 

146 


— ffH-JT'JH “ 

boxpu 

14/ 


00 58 

80XP0 

140 

' 

5g'T< k rr.’Ki^TTK »-FTr, sv r - - 

BDXPO 

149 


59 CONTINUE 

eoxpo 

150 

c 


BOXPO 

151 

c 

SOLUTION OF ECUATIONS 

SOXPO 

152 


"ETr continue 

”FtixFir~ 

■■ 15 ^' 

c 


BOXPO 

154 


TFTjRlTE.tQ.O) GO TO 72 

BOXPO 

155 


MfUie I 

eoxpo 

156 


DO '70' K*r»2 

BOXPO 

157 


Nl = l 

BOXPO 

158 

■■ 

' N^»Ni*3 - — . ... 

BOXPO 

159 


66 00 67 J'l.JE 

BOXPO 

160 


67 HRITEUH*i69) IJ,N,f fK,J,N},N = Nl,N2) 

BOXPO 

161 


IF(N2.0£,JE1 CC TO cjd 

EOXPO 

162 


hl=N2-*l 

BOXPO 

163 


N2=Nl+3 

BOXPO 

164 


■ “gc Tu do 

BOXPO 

165 


6H CCNTiNUc 

a.oxpc 

166 


CO u9 J*l,jc 

BOXPO 

167 


HKITE(rw,17Cl J,I) 

BOXPO 

168 


69 continue 

BOXPO 

169 


7C CONTINUb 

BOXPO 

170 


157 FOKNAT ( IHO.SX , 12 ,l >.1-Td C^Lo^^) 

BOXPO 

171 


168 F ORhAK IHl, lu;t ,6-K IClEv’ .-AT'.:* (A) > {A)*(X)=(8) FOR VELO 

BOXPO 

172 


iClTY POTENTIAL IX) ALONG EACH COLUMN//)' 

BOXPO 

17-3 


169 F0RMAT(8X, A(2> ,lHt , 12,1H, ,I2,2H) «E13.6)) 

BOXPO 

174 


170 FCRMATl9X,2hS t ,I2,1H, ,12 ,2H)-,E13.6) 

BOXPO 

175 


72 continue 

BOXPO 

176 

c 


BOXPO 

177 


K = NSiMEC(NB,JE,l,T,Sii, JS,1) » 

BOXPO 

178 


IFlX-NE.l) GO TO 9a 

BOXPO 

179 

c 

COMPUTE HAKE POTENTIALS — 

BOXPO 

ISO 

c 

- PHHTE)*EXP<-IK*lX-XTEn 

BOXPO 

181 


IFtX.LL.lAMl-cj) GO t6'85 ' ~ ' ' 

"BOXPO 

182 


JS = 1 

BOXPO 

183 


75 Y=YY<JS) 

BOXPO 

184 


IF tx.LT.tXTEl J5>-CH) J CO TO 82 

BOXPO 

185 


IF lX.CE.lXni JSI^uril) GO TO 82 

BOXPO 

186 


IF tX.GT.XTEl JS) ) GO Tj 76 

BOXPO 

187 


PTR-S (i , JS, I 1 

BOXPO 

is 8 


. PTI=5l2,JS,n 

BOXPO 

189 


KK = 1 

BOXPO 

190 


XB=X+D 

BOXPO 

191 


CO TO 77 

BOXPO 

192 


76 P«(X-XTE(JS) I/C 

BOXPO 

193 


PTR = ( l.Q-P ) «S 1 1, J5, I) *P *S1 1, JS, I-l) 

BOXPO 

194 


PT1*{ I.O-PIYSIZ, Ji,l )+P»5(2, JS, 1-1) 

BOXPO 

195 


KK = 0 ' , . . . . 

BOXPO 

196 


xe^x 

BOXPO 

197- 


77 CONTINUE ■ • 

BOXPO 

198 


80 IF (Xe.GT. 1.0) GO TO 82 

BOXPO 

199 


XH=(Xa-XT£ ( JS ) )*Ch 

BOXPO 

200 


IK»i=I+KA 

BOXPO 

201 


Sli,JS,IKH)=PTR9C0S{XNJ*pfl^SlNlXH) ' 

BOXPO ■ 

. 202 


S12,JS, UH) = PT I*COSlXN)-PTR«SIN(XH) 

BOXPO 

203 


KK^KiC + 1 

BOXPO 

20V 


xe=xb+o 

BOXPO 

205 


GO to 30 ‘ ' ■ ■ 

BOXPO 

. 206 


82 JS=JS+1 

BOXPO 

207 


IFiI jS.Lt. JL) GC TO 75 ' ' ' 

BOXPO 

208 


85 CONTINUE 

BOXPO 

209 


90 Il«Il+l 

b6xpo 

“irsT 
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RETURN _ __ SOXPO 

'95 vr'jtE ' ' ■ ■ '< BOXPO 

STOP _ B OXPO 

100 H?k>^A 7 Utll»ECx.‘i 2 HlhFLU£NCE COEfFICIEKT (RE At AND 1 HAG 1 HARY 1 ///I BBXPO 
110 F 0 RHAT( 1 H 0 , 3 < tl 3 « 3 dH- 6 <JX SEPARATION IN CHORCHISE DIRECTION). dOXPO 

120 FORhATdH f A( lStlP 2 El 3.61 J ' ^ ^ BOXPO 

lA O FC RH ATC 1 hQ 1 Cx« 59 HS 0 LUT 10 N OF SIH ULTAHEOUS EQUATIONS FOR THE POTENT BOXPO 
UAL FAILED! ' BOXPO 

hJIC BOXPO 


2U 

212 “ 

_213 

2f4 

215 

216 

217 

218 
219 
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•SUeSOUTINb fubCI(XX,VY,S.SfDy,3,FDY,SFDH,KSF0»KSFS,XLE,XTE»YHAX 

FOflCI 

2 


1 - • — - ,'amA‘, J.-iAX iHLT,L,NEH.H2,HB,HD.NB,HM,gK,T> . 

■■■PSRII 

3 


E IMtrtSinn nLT (2,M3) ,XLl:(rtB> .XTt I*<B) ,XX'(MdJ , YY(HB) ,S(2,MB,HBJ. 

. FORCI 

• 4 


1 .5F5X CNc ,rt 3 ) .SKJY(NB.nOJ «SFOM(fiH,HD)iK5FO(MO) »T(NM«1) 

FO'RCI 

5 


2 ■ 

FORCI 

6 



FORCI 

7 


N = 0 

FORCI 

8 


MC SHOjfriiNO Ih SUaHOUllNE SMOdTH, tSPLNl) 

FORCI 

9 


C-XX<2I-X*{ 11 

FORCI . 

io 


Hi 5*^ -J 1 *1< ) 

FORCI 

11 


I f 1 -<5FJM2) .::C.O. AMD.rlt:W.E0.2) H2=KSFS 

FORCI 

iZ 


L'L' 150 J = l,JttAK 

FORCI 

13 


HCTl=t1LTU.J), 

FORCI 

14 


N2=«Lf (2.J)-rlL'71*‘2 

FORCI ' 

' 15“ 


sC^;«2»l 

FORCI 

16 


M*l 

PORCI 

17 


K£ = <i 

FORCI 

18 


KC = 0 

FORCI '. 

19 


T{l,i)=XTE(J) 

FORCI 

20 

'.c 

■ 

■ FORCr 

“Zl' 


Til. 21=xt.tlJl 

FORCI 

22 


T( = 

FORCI 

23 


kA = ' 

FORCI 

24 


K b 

FORCI 

25 


i h ( ^l • 1 0. 3 • Cr • NZ A c U • J 00 T*J 

FORCI 

26 


iFi»:.cj'.lT aC 

" FORCI 

27 


IXCLJC^ LtACIhC fcOJi. PUINT IF IT IS TOO CLOSE TO FIRST BOX 

FORCI 

28 


!F (j).OT.0.2«D) CO to a5 

FORCI 

29 


n A *■ 1 

FORCI’ 

30 



FORCI 

31 


K C = 1 

FORCI 

32 


cosm.Lt ■ ■ ■■ ■' " 

FORCI 

33 


53 x=R.t,KE 

FORCI . 

34 


Tf' . (XM 

FORCI • 

35 


i F ( rl .t 0. 3) CC TO 30 

FORCI 

3b 


T(K,k2>=S(K.1»J,KK1' 

FORCI 

37 

so 

K K =K^ 1 

FORCI 

38 


1 f- ( 0 • t U,, ) 00 T U i 

■ FORCT ' " ■ 

39 


lF(i.l.tJ.3) at TO oO 

FORCI 

40 


CALL SPLNUl,TH,l),.Tll,&)',T(l,9),kB,T(l, 2 ) ,Tf 1 ,<2 ) »Tl 1,10),N3 

FORCI • 

41 


TINO,x2)=TU.6J 

FORCI 

42 


IFlKC.Nt.f) O'C TO' 51 

FORCI ' 

43 


kC = 2 

•FORCI 

44 


ctr TC'Tir - 

■* FTTRCT ■ 

A5 

SI 

CCNTI.Sue 

FORCI 

46 


• * “ " ” ’ 

FORCI ■' 

47 


1FIS2.0T.2) GC TO 52 

FORCI 

48 


AOJUSTHENT for TdAILIMG EllGt VELOCITY 'POTENTIAL 

FORCI 

49 


KK^rtLTt 2» J) 

FORCI 

50 


DL’N=SdRt,l fxTFl Jl-XLbl J) ) / (XXikK »-kLE< JJ H 

FORCI 

51 


T IN0,K2 J^COfAT f KZ.K2J 

FORCI 

52 

52 

CONTINUE 

FORCI ■ ‘ 

53 


* 

FORCI 

54 


kC=0 - -- 

F0RC1“ 

■^5 


lFtKl.£Q.2) GO TO 55 

FORCI 

56- 


Kl=2 

FORCr 

57 


X2 = 6 

. FORCI 

- 58 


GC TQ AO 

FORCI' . ' 

■ 59.,. 

TT 

kl = 3 

FORCI 

60 


GO TO AO 

FORCI 

61 

hO 

TIN«. 7)=xTf(J) 

. FORCI 

62 


IF IASS IXLEl J )-XX MLUn .GT.l.t-OM GC TO o7 

FORCI 

63 


btAblNG c1)GE isb FIRST b.lX CO INC Kb 

FORCI 

64 


■ NZ=nZ-1 

FORCI 

65 


hC=N2-l 

FORCI 

66 


CC o5 K=2. NC 

FORCI 

67 


TIX-, 2)=t'(K«l, 21 

FORCI 

68 


T(K^ A)=T(K*1, A1 

FORCI 

69 


T<K» cJ-Tifk + l^ 

FORCI 

70 


58 
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lf(N£w.L^.l) IC Tj 

FORCI 

71 


TlK, i;> =T t^♦x ,1:;) 

FORCI 

72 

t*> 

C <’"*T 

FORCI 

73 

r ’ 

'. -SflSUi. 

FORCI 

74 

c 


FORCI 

75 

c 

lNT;-^?LLATt. jiFi.h:iIJN aT ?C> CtNTf*^* LtAOIhG AND TKAILING EOGES 

FORCI 

76 


CALL SURF2 (ry(^) ,T( i, . J , 1 , I ( 1 , 1 3 ) . T ( 1 , 14 ) iL)U«» SFOX (i ,M2 1 

FORCI 

77 

i 

- SF JY( i,.-'2 ) .ir ,iiU ,r- 3 ) ,^ 2 ) 

FORCI 

78 

c 


FORCI 

79 


L’G 7; K-1, 

FORCI 

* 80 



FORCI 

81 

76 

TlK.lfc)^ H^,l«.)?TlAt5l-LK*T(K,U)tT(K,4) 

FORCI ■ 

82 


CALL iNToub f : , L . 1 *TI J, W) ,T lj,lbl» n 1, 19) ,2,NQ,HZ»NM) 

FORCI 

83 


CU- = T j. 

FORCI 

84 


7-( J,17) = I j?„on-n J ,17) 

FORCI 

85 


T( J ,18) =T<NC^t )=>3UH-T tJ , lb) 

FORCI 

86 

15C 

CONTINUE 

FORCI 

87 

C 

Pt7F3RT 3PA^HI5t LNlt jRAT ION 

FORCI 

88 


SC- JhAK *2 

FORCI 

89 


M-17 

TorcT 

90 


K2-=13 

FORCI 

91 


T(1,1)=0.0 

FORCI 

92 

lo6 

DC 170 K»l,Jr!AK 

FORCI 

93 

17C 

T(k, 2)=yy{K) 

FORCI 

94 


S2= JMAK 

FORCI 

95 


CALL SFCn 1(1,T(1,1) ,T(l,b) ,T(1,'/) ,NZ,T(1,2) ,T(l,iCl),Tll,lO),Nl 

FORCI 

96 


T ts j,K2)=0.fl 

FORCI 

97 


T( 1,k2)=T<1,6) 

FORCI 

98 


OC 175 K=1,JMAK 

FORCI 

99 

175 

T[X-,l ,K2)=T(K,K1) 

FORCI 

100 


iFiKi.ea.is)- cc TO lao 

FORCI 

101 


aI*18 

FORCI 

nfif 


A 2- Lo 

FORCI 

103 


GO TO 165 

FORCI 

104 

ISO 

CC^Tl^L»£ 

FORCI 

105 


KZ* JrtAK ♦! 

FORCI 

106 


T( 1,2)=0.0 

FORCI 

107 


T()*Q,2)-yf<AX 

FORCI 

Tff5“ 


OC iS5 K=i,JhAK 

FORCI 

109 

185 

T[K+1,21»YYIK) 

FORCI 

no 

c 

MING TIF COftftECTlOH 

FORCI 

lU 


K1«0 

FORCI 

112 


xCsHhqjZi-tifz-IiZ) 

FORCI 

Z13 


PR-TCNZ-1,15) 

FORCI 

'114 


PI«T(NZ-1,16). 

FORCI 

ns 


1FITINZ,2')+.5*C.LT.T(N0,2J) GO TO 195 

F^Z 


TfC 

3UHT5TJRTTTr(KC,2J-Il),Z,ZJ)ZXO) = 

TEBtel 

-ii?r 


TCHZV15I«DUNAPR 

Foaci 

..tter 


“TlMZ,16 I=DuH»PI ” 


riav 


IFIKX.EQ.l} GC TO 196 


^0. 

195 

IF tYYt4«AIC)*1.05»0.GE.YH*X) GO TO 196 ^ 

FOa^ 



Kl^=l 

F^RCr >. 

~ur 



■«Rcr^ 




Foad'- 

li<F 


TIN0,ZI-TING-1,2) 

FORCI 



TINZ,2»-T(NZ-1,2)*D 

FORCZ 

£20 


TINQ,15)«T<N0-1,15) 

TKBCr 

l^r 


TINQ, 1*J-TINQ-1,16) 

FORCI 

no 


GO TO 190 

TOTCI 

EPT 

' lit 

CCSTIhOE 

FORCI 

139 


I F1Xb5 ITTN o', 2 }-T(NZ,2)).CE.1.E-05) GO TO l9^ 

FORCI • 

'131 


NQtNO-1 

FORCI 

132 


NZ’NZ-1 

FORCI 

133“ 


TIHO, 2)»TO)0-*J., 2) 

FQBCl 

134 


■tiwo;t5 }-tt>iqtivt5 1 ~ : 

FOfiCI 

— 135" 


T<NC,16)=TtNQ.*l,16) • ■ ■ 

FORCI 

136 

,"T^ 

CONTINUE 

FORCI 

137 

C 


FORCI 

138 


CALL 1 NTGL ( T ( 1 ,2 1 Vm-VlS ) , TTI , I > ,T ( 2, 1 ) ,TTl,T91 ,2, NffiHZ ,HM ) 

TDRCl 

139 

* 

RETURN 

FORCI 

140 


m FDRCl I5T 


59 



o 


c 


c 

c 

c 


c 



SUBROUTINE MRE0C0A,T,NN,NB»KSFN»SfNK,SF«Y,5fMH»IH,'l«UNTJ 

HRED 

r 


SPLlNE-SURf ACE FIT JF MACH' NUMBER 

HREO 

3 


UlNlNSlun t}Am,I(NN«ll«SFNX(l) tSFNT f 1 ) ,SFNHf 1 ) 

MREO 

\ 


C0NST»0. 28571^29 

HRED 

5 


KSFH^DA(971 

HRED 

6 


IFfKSFh) 80,50,10 

HREO 

7 



HREO 

8 


FITTING OF GIVEN PRESSURE/NACH TO A SFLIHE-SURFACE 

HREO 

.9 



HREO 

10 

10 

IFINB-KSFM) 80,15,15 

HREO 

11 

15 

CONTINUE 

HREO 

12 


KP»701 

HREO 

13 


DC 30 IP=1,KSFP 

HRED 

19 


SFHX( IP I-OACKP )/DA<2A) 

HREO 

15 


SFNVIIPl=DAfKP*l»/0A(2A) 

HRED 

16 


SFHHUP )=DAIKP*4 ) 

HREO 

17 


DAC96)=1, INPUT DATA ARE PRESSURE COEFFICIENT 

HREO 

18 


DA{96>=2, INPUT DATA ARE LOCAL MACH NUMBER 

HREO 

19 


IF«OA(96)-1.0) 75,20,25 

HREO 

20 


CONVERT PPESSLkE COEFFICIENT INTO LOCAL HACH NUMBER 

HRED 

21 

20 

SFMH< 1P) = SCRT(5.A( I.2/J ( 1 . +0. 7*SFMHI IP ) )*«CCNST 1-1. 11 

HRED 

22 

25 

continue 

HREO 

23 

30 

KP=nP*3 

HREO 

24 

. 

SPLINC-SURFACE FITTING OF DATA 

HREO 

25 

•1^0 

CONTINUE 

HREO 

26 


IF t IPRINT.NE.C ) WRITE! IH, 1001 

HREO 

27 


CALL SJRFKNP ,KSFN,T,SFNX,SFMY, SFMH,IPRINT1 

HREO 

28 


Rt TURN 

HREO 

29 



HREO 

30 


presently INPLt CF pressure COEFFICIENT IN 

MREO 

31 


A POLYNOMAL FORM IS NOT ALLOWED 

HREO 

32 


THE FCLLOKING IS FUF NACH INPUT AS A POLYNOMIAL M=A0-»A1« 

HRED 

“35~ 

50 

CCNT INOE 

HREO 

34 


NSf .•'=0 

TSRED 

35 


SFM»( 1)=JA( 71) 

HREO 

36 


SFwH(2) =DA( 72)RDA(2A) 

HREO 

37 


SFHHI J) =C.O 

HRED 

38 


vC To ^0 

HRED 

5^ 



MREO 

40 

75 

IPR=9b 

HRED ‘ 

' '4r 


GO 70 85 

HRED 

42 

3C 

IPR=97 

HREO 

'43 

S5 

WRITE 

HREO 

44 

C fIrt.liOllt’K 

HREO” 

^ 


STCP 

MREO 

46 

IOC 

FORNAT(1HO,10X,73HCONPUTED NACH(X,Y) - A0+A1*X+A2*Y* SUM OF HID* 

HREO 

47“ 

t(R(U*«2)*tALCG<RtI)**2)n 

HREO 

48 

110 

FORWAT(lhO,10X,l<iHNRED BAD DATA, 15) 

HRED 

49 



ENO 

HRED 

50 
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^ QcrAtm 




SUBROOTIRE INTGLi X«T«VR*V|«S«N«^0«ttZ«^) ' ' 


c 


IMTEGRATIOM BASED OH SPLINE FUNCTION -i.. 





OINEHSlOH Xlli,r<NN,2»,SCNN,ll ■ • ' 



c 


DEFINE HJ) 


• -# 



DO 20 1=2, NO 


» 


20 

S(I*1I=XI1)-XII-1) 


* 



IF INQ.EQ.2> GO TO 50 

- 

• 

c 


DEFINE TRI-OIACONAL COEFFICIENT MATRIX' 


V 



DC 25 I«2»NZ 

IttfCL . 

X9 



S(I,2>= S(I,I)/6.0 

INT«L 

ti 



S{l,3)=(SCI,l)*S{i»l,l])/3;o 

INTes. 

12 


25 

SII,4)= S(I'<1, 11/6.0 

IKTCi 

13 



SI 2,2)=0.0 

IHTCL 

14 



SINZ, 41=0.0 

INTO. 

IS 

c 


DEFINE RIGHT-HANO-SIOE COLUHM MATRIX 

IMTa. 

.16 



K-1 

INTCL 



35 

DO 40 1=2, MO 

IHTGL 

IS 


^0 

Sn,5) = ITU,K)-Y(I-l,K)l/S<I,l) 

IKreL 

IS 



DC 45 1=2, NZ 

I«T6L 

2S 


45 

SII,6l> $n-»i,si-su,5J 

“THTH 

21 

c 


SOLVE FOR COEFFICIENTS OF SPLINE FUMCTION M(4) 

IMTSL 

22 



CALL tRlbl(2,f>Z,S(l,2).SIl,3),S(i,4),Sll,6l*S(l,7),Sll,8l,S(l*9ll 

IMT^ 

2S 


50 

CONTINUE 

INTGt. 

24 



S( 1,71=0.0 

INTCL ^ 

25 



S (NO, 71=0.0 

IMTCt 

26 



VI»0.0 

■“ftTRE 

if 



OC 60 I =2, NO 

IMTCL 

28 



VI-V1*0.5*SU ,11 

INTCL 

24 



i ♦( Y(1,K1 *7(1-1, Kl-S( I, 11 *5(1 ,11* (SI 1,71 ♦SI 1-1,711/12.1 

INTCi. 

30 


60 

CONTINUE 

INTCL 

31 



IFIK.EQ.21 RETURN 

INTCL 

32 



VR=VI 

IMTO. 

33 



IFIN.EO.il RETURN 

INTCL 

34 



R=2 

nntc 

35 



CO TO 35 

INTCL 

36 



END 

INTCL 

3? 


3UBi<QUT!NE Tk IC 1 1 KI*K3« A , B«C«D» V ,£ >F > 

01 hfc'«SlT3S“Ttl J7Bf TlTC'a J 7DI n , V U1 ,t < nVf VIT 

IF (K3.N E.K1 ) GO TO 5 

V(K3") = 0'fKT»7B(K3T 
Tcc TUk N 

5 CCnTINOr ■ ■ 

E(KI)=B(k1) 

■'f ( K l)-U<MT/t(Kil 
H^=^l+l 

cr 11 l=«.3,Kj 

t { 1 1 = 3(1 )-A( n*C(I-l) /£U-1) 

1 C F(i)=(un)-Atn*F(i-r))/€m 
V (K jl=F t\ 3 J 
it r =Tn -TT 
L'C 3C J = l,^2 
I=K 3 -j “ 

20 v< I 1 =F ( I )-C n 1 »VU_*U/E<J,I _ _ 

FCT'J^K ' ■■ ■ - — - - 

fc Nt» 


- ^{01 2 

■TRIM ar 

TRIOI 4 

TRiOl W 

TRIOI 6 


TRIOI 8 

nrroi 9~ 

TRIOI 10 

TfflCI 11 

TRIOI 12 

TRTXn 13 

TRIOI 14 

"■ TftlDl IT- 

TRIOI 16 

TRitfr 17 

TRIOI 18 

” rmn i9 

TRIOI 20 
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'Onlo* rjrj r> o o o o o 


c 



SUBROUTINE SPLNUHI , X, Y,DY,N,XX , YY,DYY<NS«OSI. 



2 


X , Y'iTTYilNTfc kPCLAt JON INUfcPtNOEHr, DEPENDENT VARIAbLbS*. AHU 

OV/Ux 




Ni =No. OF Interpolation points' 



4 


XX, yy =ispuT Ino^'pendent and' dSpendent variables 

* ~ 

SP'CMl 

“5 


N -=no. of Input points 


'SPLNl 

6 


OYY =D«YY)/CIXX) FOR' INPUT'DATA 


SPCHT ■ 

“7 

. 

NSHOs =c-.MfC'..s ;’F smu-tthing and pre-interpolation 


SPLNl 

-8 


Ilr.NSioN X 1 I ) ,Y( ; 1 ) ,XM1 ) .'YY ( n iXiYYIl) ' 


5TOST“ 



Ni‘ = i> 


SPLMl 

10 


I P'-.^ =N5«3S/10 

' 

StfUBT"' 

■ n 


nS“L=nsH0S-IPPE»10 


SPLHl 

12 


IFl IPRc.tC.O) GO Tw 3J ' ' ' 


SPLNl 

13 


STORE INPUT jiTA FJR P RE - 1 Nl£k P OL * T I ON 


SPLNl 

14 


Ij C 2*^ J ■= 1 , N 


-SPtMl" 

■T5" 


Y{ J) =XaIJ) 


SPLNl 

16 

20 

DY(Jl*YYtj) 


SPLHl 

17 


CALI .SPI SET IN, XX, YY.DYY, 0.0,0) 


SPLHl 

18 


PRI-INTERPULATIJn 


SPLHl ■ 

19 


DO 25 J=2,N 


SPLNl 

20 


ri*T^-7 ' 


• spcra ■- 

71 


12=1 1>1 


SPLNl 

22 


sxU2) = Y(i> 


SPLNl 

23 


YY(I2)=L.Ym 


SPLNl 

24 


X> ( ! 1 ) = 0. 5A (V < ( II.) 


SPLNl 

26 

25 

CAUi. SFi.N2(XA<Il),i,U,Y,0Y ,0YY, YYMU , 0UM,1 1 


SPLNl 

26 


N2=7-N-l 


SPLNl 

27 

30 

CGNT INUu 


SPLNl 

28 


IE INSFD.EL.C) GO TO 50 


SPLNl 

29 


SFOCTF Input CATA XA.i YY 


SPLNl 

30 


C'AL^ SnCOTmI Nl , RX , Y.Y,r ,NS(*0) 


SPLNl 

31 


iNTErP'lLATE Y AT / F kO -i XX, YY, 3YY AND CALCULATE 0Y = 0(Y)/DIX) 

SPLNl 

32 

50, 

CALL SPls7TCN2',X5tT77,757T',0.T5T01 ' 

" — — — 

“SPLHl 

33 


CALL SPLN2 (X ,M ,n2,XX, YY,DYY,Y, DY,2 ) 


SPLHl 

34 


RETURN 


SPLNl 

35 


LnC 


SPLNl 

36 


SUaKG UTINb SPl.h3(XP,NP,N>Xty^0»SPf «SPP>K) 

■'DlhENSION X<ll ,Y(U,0(l>«XPm,SPFm,SPDIl) ■ ■ ~ 
EVALU ATES A hjlTURAL C UBIC SPLINE AMOX'^S FIR ST OE RIVATI 
SLOPE ARRAY D CALCULAtSD BY SPIs“e'T ANO USING THE INPUT* 
AkRAYS X AND Y 
DO 

GO TO b 


VE USING 
DATA 


10 J=1»NF 
I F(XPU).LT.X(1I.0R.N . £Q,H 


00 2 1=2, N 

1F(XPU>.LT-X(ln GO TQ_^ 
CONTINUE 

SPF (J) = Y(N)0(X>_?l^>‘‘i,lZ?<tN>> 
IFIK.EQ.II Go'TU lo 
SPDI J )=D(N) 


5PLN2 

SPLN2' 

5Pyi2 
SpLh2‘ 
SM-*« 
SPLN2 ■ 
SPLN2 


2 

3 

5 

7 

6 


GO TO 10 

Cl = l-/(X( I J-X^I-X) J 
LA=*i': )-xp(j) 
C3=xPUI-X(I-l ) 
C^=C2*C1 
C5=C3*Ci 


SPF (Jl = Cb«C3*l (l. + 2.«CA)«Yn )-C2<=DI I J ) 
i *C5J*yU-ll«C3*DI i-lJ> 

iF(K.FC.ij go' fd icT “ 

C6=2. AC2-C3 

C7=2 . AC'3-C2‘ ■ ■ - - - 

SPS(J> = C1«C1»(C3»(2«*(1.*C1»C6> »Yn>-Cfe*t)II)l 


10 


I -C 2 * ( 2 .*U. ♦C 1 *C 7 J*YI I - 1 ) + C 7 *D I I-l» H 

GO TO 10^ • 

SPF fjy=Yfl)-u (XK'fxI 1)-XP< jn ■ " 

IF(K.tC.l)_ uu TO 10 

5PC(jr=eu)' - - 

CCNTI'-L. 

RETUVN ~ ■ 

cNC 


SPLN2 

9 

SPLN2 

10 

SPLN2 

11 

SPLM2 

12 

SPLN2 

13 

SPLH2 

14 

SPLN2 

15 

SPLH2 

16 

SPLN2 

'17 

SPLN2 

18 

SPLH2 

19 

SPLH2 

20 

SPLN2 

21 

SPLN2 

22 

SPLH2 

■ 

SPLH2 

.-24>V 

SPLH2 

25 ■ 

SPLN2 

26 

SPLN2 

rr 

SPLH2 

28 

SPLN2 

- “ '2'9 

SPLN2 

30 

Spln2“ 

31 

SPLN2 

32 

TPtH2‘ 

33 

SPLN2 

34 
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ORIGINAL BAG£; 

OF POOR quality. 



SUS»OUTIhE SPlSETiN,7,y,:,KMS, ifr«sj 

SPISET 

2 


X ( 1 ) 

SPISET 

3 


DATA JMA/itti>‘/fcrfi*io33r;.o7»l«j71?9677/ 

SPISET 

•<( 


om-=^o. 

SPISET 

5 


IF(N.cO.l) i>lTUAn 

SPISET 

6 


DCN)=0. 

SPISET 

7- 


an=n 

SPISET 

8 


DLj i 

SPISET 

9 


I=h*2-II 

SPISET 

10 



SPISET 

11 

i 

Y( ; ) = ( y ( i)-Y{ 1 - i ) j // ( ; ) 

SPISET 

12 


ll-(s.cC.2) Cj 7., ; 

SPISET 

13 


OC 2 II=3tN 

SPISET 

14 


TiK7 J-I f 

SPISET 

15 


i3( Y 1 1 ]-y( i-ll ) / (XU )*x 1 1-l) ) 

SPISET 

16 


Y n ) = l.S<=C f !-l 1 

SPISET 

17 

2 


SPISET 

18 


1,2 3 J-l>J9ix 

SPISET 

19 


LC 3 

SPISET 

20 

3 

;i!-r5=HA(Y(i)->(i (i-2)-(D.3-x(i))«Dn)i-<K-i.o»#on-ij 

SPISET 

21 


Jl' A !=j.N 

SPISET 

22 


X a ( : -1 ) c .‘/X 1 i )-*. ) 

SPISET 

23 


Y t : I - . . j JJi J ; j -ay ( i J A( < ( ; ) +A( 1-1) ) tYII-l ) 

SPISET 

24 


= 12) ) 

SPISET 

25 


3C ; : = 

SPISET 

26 


( n 

SPISET 

27 


‘ji»Tt:)-;v(, .A.ii-.)».'(i)) 

SPISET 

28 


(i-i )-»2.*o( 1 1) 

SPISET 

29 


_tl-l )='u.b*l3A’f; ♦xl) 

SPISET 

30 

6 

SaV : ^ l3 

SPISET 

31 


L ( '^ ) X S »,V t. 

SPISET 

32 


IF(7r^'i5.A' .11 -33 T7 3 

SPISET 

33 


A 7S =.. . 

SPISET 

34 


If (N.tC.2 1 CD TO ~ 

SPISET 

35 


jfi 7 1=3, N 

SPISET 

36 


C=2.ACII-1) 

SPISET 

37 

7 

P)<S = »lMS+t tC-YDl I-2)-3.7Y n-i) ) / X ( I-l)t (D( I J+C-3.»Y( I ) )/X<I ))**2 

SPISET 

38 


TfHr=rr?Ti5Fn uTi/Ah i ----- - - 

“ SPISET -- 

“3T 

3 

00 9 1=2, H 

SPISET 

40 


Y 1 1 ) = Y( i)*x n )+y n-u 

SPISET 

41 

9 

Xt E1 = XUJ*>U-1) 

SPISET 

42 


. Tua N 

SPISET' 

43 


END 

SPISET 

44 
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SUBROUTIHt SMOOTH IH.X»Y,T»MSHOSl 

tSNtSld 


5 

THE Y ARRAY IS SKOOTHOO BY A LOC«- Fl« 

CUBIC HEIGHTEO BY M 

rOlNT LEWT SOARES 

SKSOTM 


■ 

DKSNSION XfX}<YUt«T(.U 
lFfH.i.T.S> RETURN 


mmm 

B 

DO 10 «S«1»MSM0$ 
TI11>NS 



AN-M 

S»<TIH*CXIMI-X<1J I/AN1**2 







DO 4 L>1«N 

K>HIH0tH>4«HAXQ(liL~2)l 




DO 1 I°l*20 

: 



1 TJ1»*0. 

DO 3 M«K«K4 




Wl./»S*IXIU-XtH) »**2I 
R«1.0 

smssur 



00 3 J-l»4 
IA-I-4 

StISUTH 

USSQtH 

2 * 


RR-1-0 
00 2 J>1«4 

SmJSiW 

SftQOTH 

ZZ 


J4-4*J*14 

T< J4J-T(J4>*R4RR»R 

' SMOOTH 

SMOOTH 

23 . 

24 


2 RR-RR*X<H) 

TtI+16l-TiI*16)+R*YCMI*»l 

SKOQTH 

SMOOTH 

25 

26 


3 R*R«XIH> 

CALL CHLSKYJT»4,T117>»1»41 

SMOOTH 

27 

26 


M«L-UL-1J/5I45 
IFCL.GT.51 T(L-51»T1M*20) 

SMOOTH 

SMOOTH 

29 

30 


- TtM*20»-0. 

R-1.0 

SMOOTH 

SMOOTH 

31 

32 


DO 4 J=l,4 

T(«*2CI=TtH+2C HR*TJ J*16» ' 

SMOOTH 

SMOOTH 

33 

34 

“ 

36 

37 

38 


4 R-R4XCL) 
I4=N-S 

SMOOTH 

SMOOTH 


DO 5 L=i»5 
HL=M*L-( 

SMOOTH 

SMOOTH 


J4*=I4*L 

5 Yt J4)=*TtML + 20 1 

SMOOTH 

SMOOTH 

' 39 

40 

41 

42 
43' 


10 CONTINUE 
RETURN 
END ■" 

smooth 

SMOOTH 

■" SMOOTH 



_ . ^ 
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r> 


ORIGINAL BAG'S IS 
OF POOR aUALITV 


SUBROUTINE CHLSKYiA.HVB .HfNX) 

DIMENSION AIK)(,1» iBINX,!) ' - -- - 

CHOLESKY DECOhPUSITlON IS USED TO SOLVE THE MATRIX EQUATION AX»8 
MHEKE THE CCEFFIOIENT MATRIX, A, IS SYMMETRIC. OM OUTPUT" X IS 
STORED IH B 

IFIN.EQ.II GO 7b"6 - . - - - - 

DO 2 1=2, N 


CHLSKY 

CHLSKY' 

CHLSKY 

CHLSKY 

CHLSKY 

CHLSKY 

CHLSKY 


Il-I-l 
00 2 J=Ii 


CHLSKY 

CHLSKY 


♦ 

5 

« 

7 
8_ 
— q 
10 



OC 2 L=1,I1 

CHLSKY 

11 

2 

Att,J) = A(I,J)-A(L,lI*ML,Jl/A(L,L) 

CHLSKY 

12 


00 5 K=1,M 

CfCSKY 

13 


00 3 1=2, N 

CHLSKY 

lA 


I1=I~I 

CHCSKY 

TT' 


DO 3 L=1,I1 

CHLSKY 

ifi 

.. . 

B(I,ltr=glf,Ki-A(L*'i )*STCVK)/ACL,L) - 

~CTORY 

- X7 


DO A 1=2, N 

CHLSKY 

10 

“ — - — 

”11^-1 --- 

cra-iKY 

'19 


00 A L=l,Ii 

CHLSKY 

20 


KI=H-'I1“ ■ 

CHLSKY 

"21 


NL=N+ 1-L 

CHLSKY 

22 

' A 

6 ffTHK ! =TTh I , K )-A ( Nl , NLl *e ( NL , >TJ7miL , HL 1 

CHLSKY 

23 


OC 5 1=1, N 

CHLSKY 

2A 

5 

B('i,Kj=R'r,K)/ATi ,n ■■ ' - -- 

■ CHCSKY 

25 


RETURN 

CHLSKY 

26 

6 

Ati,i)=r;/Ari,i) 

CHLSKY 

Z1 


CD 7 L=1»M 

CHLSKY 

28 

7 

s(i';cT=iTr,T.TY?n,L"T — - 

— CHLSKY 

” ^<3 


RETURN 

CHLSKY - 

30 
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SUBP3 JT INI SUPFlCNM.Ni T ,‘ABX,ABY fABH«IRITEl 

SUf^l 

2 


"FfT =)ATA (h PCINTi) BROUGHT THROUGH ABX,ABY,AB« 

SURFl 

3, 


I __ 1..HFCRAKY ARnAY FOR SPL I M ■■ -SURF ACE FITTING 

SURFl 

A 


; — hTY arrays for INPUT POINTS TO BE FITTED 

SORFl 

5 


- I-^/tPENOENT variable X 

SURFl 

6 


ASYtlJ - INDEPENDENT VARIABLE Y 

SURFl 

7 


Abhin - CbNES IN AS dependent variable of X AND Y 

SURFl 

8 


goes out as coefficients of spline-surface 

SURFl 

9 


LIF.'^SI^N ;<'i|-tl),A3X(l)'»ABY(l)»AeH(ll 

SURFl 

, 10 



SURFl 

11 



SURFl 

12 


NP1=N *1 

SURFl 

13 


NP2=NPif L 

SURFl 

lA 


NP3=SP2+1 

SURFl 

15 


nPA=NP3*1 

SURFl 

16 


IF (n.EO.O) gD to U ' 

SURFl 

17 


i»N , 

SURFl 

18 


T t I ,NP2 ) = AtX< I > 

SURFl 

19 


T ( I ,nP3 > = AeYl I ) 

SURFl 

20 

2 

T(:tHFA) = AbP(I 1 

SURFl" 

Z1 


-C 6 i=l,N 

SURFl 

22 


T<l,I)=0. 

SURFl 

23 


T ( : ,NPl ) = i,C 

SURFl 

2A 


T 1 NP l9i)~L«G 

SURFl 

25 


T(NP.',1 ) = T(I,FP^J 

SURFl 

26 


T INP 3, 1 )=T ( 1 ,FP3J 

SORFl 

27 


6Hi iM-i 

SORFl 

26 


Cw 1 = 1 

SURFl 

29 


iP 1=1 •» 1 

SURFl 

30 


„c t, j=;pi,N 

SURFl 

31 


XX = T t I *NP2)-T ( J,NPi ) 

SURFl 

32 


YY=T( i ,NP3)-T ( J,N?3 1 

SURFl 

33 


H = )i yfl'>A + 'rY- 

SORFl 

3A 


T< I ,J)=HiALCG(F) 

SORFl 

35 

b 

TI Jtl )=T(I .J) 

SURFl 

36 


BC 3 1=1,3 

SURFl 

37 


IPF=I ♦N 

SURFl 

38 


DO d J=1,A 

SURFl 

39 


JPN= J*N 

SURFl 

AO 

a 

T { IPS, JPn)=C. 

SURFl 

A1 


K=«5I N=B( Nrt,NF2,l,T (1,1I,T(1,NP4»1 

SURFl 

A2 


IF<K.tG.l» GO TO 9 

SURFl 

A3 


HFITE 

SURFl 

AA 


S T C P - . 

SOREl 

45 

9 

CONTINUE 

SURFl 

A6 


5TC»E JNt3)-CCeFFICI€NT IN ARRAY ’ ABH ' ' 

SORFl 

A7 


OL 12 I=>1,NF3 

SURFl 

AS 

12 

AeHm = T(I,NPA) 

SORFl 

A9 

13 

IF(UITE) 18,18,14 

SURFl 

50 

l<i 

WRITE ( IK, 200 

SoRFT" 

5T 


WRIT£CIH,12CJ IABH(U,I = NP1,NP3I 

SURFl 

52 


IFIN.EC.Ol CO TO IS 

SURFl 

S3 


DO i<> 1 = 1, N 

SURFl 

5A 

16 

HRITC I1H,110) I, ABh( I) ,AdXU I.ABYU) 

SURFl 

55 

18 

CONTINUE 

SURFl 

56 


RETURN 

SURFl 

57 

111) 

FCK*1AT(10x,I6,lP3ei4.7) 

SORFl 

58 

120 

FCRNATl62X,iP3£i4. 2( 

SURFl 

59 

2ire 

FCF*tAT rii^uTJCXiTTinrliP.l r ( 1)5 A'2 = IX-X< i ) )**2*( Y-YI I »)A*2/ 

SORFl 

60 

1 ihO.lC/ .8Ahl I^ 'JI«tHsIONLESS CL.CJKLIN4TbS - OISTfthCE/ChORO LENGTh)/ 

SURFl 

61 

21hO» 13* 1 If I I)fl0X»4>^>(I),10x,4HY(n» 

SURFl 

62 

315> •2^4'. , 12> tat-AltliiX 

SURFl 

• 63 

220 

F CF H A T < ji < J 2 cHC 1 0 I ’’i 

SURFl 

6A 


ENC 

SURFl 
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ORIGINAL BAGE 
OP. BOOR OJJALlTy. 


^UBPJ UTIWE SUR>-.;(£I»Z 2 tJlTJ 2 »HAif,VAU)«VLUlC»¥HWf«XI»ri.tHI«HcKI 
C " COmFuTE VAHiE OF SPHHE-5U«FACE >1TT€0 OA^ AT A 
^ Z1»Z2~ CDORiilKflrES OF THE .FOIXT WgH£ T>g 

C * ■ VM.U - FITTED VALUE SOUGHT ' 

C VLUX - GRADIENT OF FITTED OALUE IH X 

C ‘ VLUy - GRADIEM OF FITTED VALUE I« Y 

S xi«yr»Hi - ARRAYS FOR KHOHW pBOpggTgES iw sgSBZs^si^^eyoiffl’ ' '' 

C N — NUHBEK OF POINTS IN Xl« VI AfiftAVS - 

C MXY«0 X=-iHJIt Y»2ZUI WHERE 



NXY>1 X^ZIID* Y°Z2(J) WHERE 
NXY>2 Y>Zim« X-Z2IJI WHERE 




IP 


OlHEMSION XXflltYmt*HIIl) 

OlHENSION Zlll).Z2(Il>VALUfll«VLUXIl)9VLUTfA} 




fM 


NPl-N *1 
NP2=NP1*1 



• ■ 



NP3“MP2*1 

IFtHXY.EQ.il X>Z1(1) 



StK^ 

StH^? 

18^ 


lFtHXY.EO.2) ¥-Zltl) 


- 

svaFi 



00 40 J^J1«J2 



StaF2 

20 


IFtNXY.EQ.l) Y-Z2IJJ 




21 


1F(HXY.E0.2) X>Z2(J1 




22 


IFtHXY.NE.O) CO TO .10 




"21* 


X>Z1(J) 



^RFZ 

24 


Y»Z2(J» 



SURF2 

~Z97 

10 

CONTINUE 



SURF 2 

26 


IF IK-2) 13>12fll 



SUSF2 

Zl~ 

u 

VLUYlJ»«HItNP3> 



SURF2 

28 

12 

VLUXI Jl=Hm«P21 



SURF2 

•29 

13 

VALUIJ)*HIO>Pl)*Hl(NP214X*HI(NP3)4Y 



SURF2 



IFIN.EO.O) GO TO 40 



SURF2 

31 


DC 30 I=1*N 



SURF2 

32 


^URF 2 31 “ 

SURF2 '34 

TS 

SURF2 _ 36 

SURF2 ■ ■ 37 

SURF2 38 


21 

22 

23 

30 

40 


IF IK-2) 23,22,21 


SURF2 

- 3 ^- 

VLUYI JJ^VLUY UWHb^TY 


SURF2 

40 

VLUXI Jl-VLUXI J1+HB4TX , 


SURF2‘ 

4i 

VALUl J)=VALIjI il+HI 1 I i *H*HA 


SURF2 

42 

CONTINUE 


SURF2 

43 

CONTINUE 


SURF 2 

44 

RETURN 


SURF2 

■ 45 

End 


SURF2 

46 


TX«X-XI III 

TY-Y-YKU 

H * TX4TX*TY#fY 

_ _ 

IF(H.Gf.Ol') HAsALQGIH)" 
HB” 2 . 4 (l.+hA)*HUn ■ 
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FUNCTION CINOl«S) 

CIN 

2 

TT 



cm 

3 

C 



ClN 

4 

C 


UMir 1 f¥y>^>rfcTi r t 

”tIM 

y* 

c 



CIN 

6 

nr 



mm ^ THI 

~ 7 

c 



cm 

■« 





'■ # 



X>Xl 

cm 

to 



IF {XI 1«2«2 ^ 

^V4friH 

11 


1 

ses»se 

CIN 

12 



X<— X 

r^i 

13 


2 

X2-X*x 

cm 

14 

, 



cm 

If* 

: c 



cm 

16 

' C 



CSN 

rr 

, c 



cm 

IS 


3 

V«itfX2/4S.0-0.6i*.0$«X24^1.0)«X2/ia.O-l.Ol*X-»1.570T9633 

CIS 




U>f 1X2/45. 0»l,.0)«X2/24. 0>1.0l*X2y4.O-.57721566»-AL0£f XI 

cm 

20 




cm 

Zl 

c 



cm 

22 



«rniw:^-TOM?[:7^ii:i!»i;n:»wT3j^>F4 nr% iin:tWTM;n.4 i ^ 

CIN 

zr 

c 



CIN 

■ 24 

i 

4 

P«U iX2-»14^ 394119 )»X2«47. 411538 l«X2«8.4933^6)/( f((k2«2i. 3610551 

CIH 




1 *X2*70.376496)4X2+30.038227»*X1 

CIN 

26 










ORIGINAL BAG'S 
OP POOR aUiLITY 



FUNC-tlUN K51n£rilH>N,LiA.,B) 


.•■J 4 


OlNtNSlON. AI«.1)»8(M,11 
SC 30 I - 1-,N 




C * f . 0 

NSIHER 

5 

' 

SO 10 J = 1,H 

■ NSIHER 

6 

10 

C » AflAXl <C«AES<A [1 « J> n 

»S»SR 

?. 


IF(C.tQ.OiO) CO TO 1000 

RSltSR 

8 



H^I«R 

4 

?0 

A(I ,J) = A( I,J J/C 

nSIHER . 

10 


L.C 3 0 J = 1 1 L 

HSINER 

11 

30 

6(l,J)=BtI,Jl/0 

MSINER 

12 


IF (N.lS. i) CO TU 205 

NSlt^R 

13 


iNr = h - 1 

HSIHER 

14 


00 200 J ^ l,^^ 

NSIKER 

15 


c = 0.0 

nSIHER 

16 


^ = 0 

HSIHER 

17 


sc AO i = J»N 

HSIHER 

18 


0 =■ A3S(A(I«Jn 

HSIHER 

19 


IF (C.CE.Dl GC TO AO 

HSIHER 

20 


K » I 

(fSlHER 

“21T 


C = 0 

HSIHER 

22 

‘tO 

CCNTINUE 

HSIHER 

23 


IF(f>.EG.0.UK.C.LT.i.c-7» GO TO 1000 

HSIHER 

24 


IMK.EC.J) GC TO 70 

HSIHER 

25 


00 50 JJ = J,F 

HSIHER 

26 


C = A (J» J J } 

HSIHER 

27 



HSIHER 

28 

50 

AIK.JJ) « C 

HSIHER 

29 


-L t-.0 Jj = 1,1 

HSIHER 

30 


C = = ( J, J J ) 

HSIHER 

31 


dIJ.JJ) = B(K,JJ) 

HSIHER 

32 

fc 0 

3 <K .JJ J = C 

HSIHER 

33 

7C 

C = Al J, JJ 

HSIHER 

34 


J? = J + 1 

HSIHER 

35 


CL bO JJ = JP,^ 

HSIHER 

36 

80 

A( J, J J) = A( J , j J )/C 

HSIHER 

37 

SO 

OC IOC JJ = 1,L 

HSIHER 

38 

100 

6( J 1 JJ) - B( Jt JJ J/C 

HSIHER 

39 


00 20 C 1 = i,N 

HSIHER 

40 


IF(I.IO.J) GO TO 200 

HSIHER 

41 


C “ Ad ,J) 

' HSIHER 

42 


DC 110 JJ- » JP,N 

HSIHER 

43 

110 

A(I,JJ) = A<I,JJJ - CAA(J,JJJ 

HSIHER 

44 


OC 120 JJ = l.,L 

■HSTWES 

45 

120 
■ 200' 

cti.JJ) = 8(1, JJ) - CA3(J,JJ) 
CChfihUE 

HSIHER 

’ "hsTher' ■■ 

46 

TT 

205 

C = A(N,N) 

HSIHER 

48 


IF (A8 SIO.lt. l.E-71 GO TD lOOO 

hsIhFr" 

49 


DO 21C J = l,L 

HSIHER 

50 

210 

8(N,J) > B(N,J)/C 

HSIHER 

51 


lF(r<.EQ.l> CO TO 230 

HSIHER 

52 


DC 220 I » 1,FB 

hsih£r 

53 


C - Atl.N) • 

HSIHER 

54 


DO 220 JJ = 1,1 

HSIHER 

55 

220 

E(I,JJ) = 8(1, JJ) - CABIN, JJ). 

NSIHER 

56 

230 

ftSIHtK = 1 

HSIHER 

“ 57 “ 


SeTURH 

HSIHER 

58 

1000 

NSINER - 2 

' HSIHER 

59 


KtTUSN 

■ ■ HSTHEir 

60 


END 

HSIHER 

61 
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fUNCTION HSiHECIM»«*L*As».^.^ 

, *’ *x. ^ . 



DIISNSION A<K<l)«8fH*U 
COMPLEX A,B,C ' -r ' ' 




DO 30 I > liM - 

C • 0.0 ■ f . ■ 





10 C-*M«XHC,AB5<aEW.tA<l< 




If fC.EO. 
00 20 J 


0.01 GO TO 1000 
- 1,N 


“2o 


AfIsJI « 
DO 30 J 


Aft>J>/C 
* 1 «L 


30 


SflvJI - Bti*41/C 
Xf:tN.E««ll 60 TO 205 




W 1 » H - 1 
80 200 A « ivKIi 




C 

K 


0.0 

0 


00 40 I > J*M 
8-ABStKEAt.<Atl<JIII*A8StAIM66<»tI«JHI 



IFIC.GE.Ol GO TO 40 
K - 1 


a 


C « 0 


23 

40 

CONTIMUE 

KSIMI^ 

24 


IF<K.EQ.0.0R,C.LT.1.E-71 GO TO 1000 

KSIM^ 

25 


IF(X.eO.J) CO TO 70 

NSIHEC 

26 


00 50 JJ ' J*M 

nsifsc 

Zf 


C > AfJ.JJl 

NSIHEC 

28 


AtJfJJ} • A(K*JJI 

NSIt^C 

29 

50 

AIK.JJ} > G 

NSIHEC 

30 


00 60 JJ > 1«L 

HSIIKC 

31 


G . B(J«JJ) 

NSINEC 

32 


« BtK.JJ) 

hSTnEC 

33 

60 

BIK.JJ} » G 

NSIHEC 

34 

70 

G - 1.0/AU.Jj 

NSlNEir 

35 


JP = J * 1 

NSIHEC 

36 

- 

D^Q" 80” JJ =" JP.H 

' NSIHEC 

37 

80 

AtJ.JJI - AIJ.JJl^G 

NSIHEC 

38 

90 

00 100 JJ = 1«L 

NSIHEC 

39 

100 

BfJ.JJ) BCJ»JJ1»G 

NSIHEC 

40 


DO 200 I = l.F 

NSIHEC 

41 


IFU.EO.J) GO TO 200 

NSIHEC 

42 


G = AU » J) 

NSIHEC 

43 


DO 110 JJ = JP»N 

NSINEC 

44 

“■ ilo 

*«I»JJ) = AUjJJl -G»A<J»Jj) 

NS I NEC 

45 


OQ 120 JJ = i.L 

NSIHEC 

46 

120’ 

b(I,JJ)'= eU.jjJ - G»atJtJJ) 

NSIHEC 

47 

200 

CONTINUE 

NSIHEC 

48 

205 

G = AIHtN) 

HSIMEC 

49 


IF lASSIkcAi. IG ) ) * AsStAlNAGlG) ) .LT.l.E-7) GO TO 1000 

NSINEC 

50 

■ “ 

DO 210 J = iiL 

NSINEC 

51 

210 

E(N,J) = B(n«J)/G 

NSINEC 

52 


'IFIN.EC.D GO TO 230 

NSINEC 

53 


00 220 I = l.M* 

NSIHEC 

54 


CC 220 JJ = l,L 

NSINEC 

55 

220 

5(I.JJ) = 3ti*jJ) - AH,N)«olN,Jj) 

NSINEC 

56 

230 

HSIMEC = 1 

NSINEC 

57 


UtTOsN 

NSINEC 

58 

1000 

NSIHEC = 2 

NSINEC 

59 




60 


ENO 

HSIMEC 
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